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SOME ELEMENTARY PROPERTIES OF ILL CONDITIONED 
MATRICES AND LINEAR EQUATIONS 


N. S. MENDELSOHN, University of Manitoba and Canadian Armament 
Research and Development Establishment 


Introduction. It is a well known phenomenon that for certain sets of linear 
equations that errors due to rounding off during the computation can lead to large 
errors in the final solution. There have recently appeared in the literature a 
number of papers devoted to this phenomenon, usually called “ill condition”, 
mostly concerned with methods of keeping computational errors in control. The 
devices used are either special ways of solving the equations or of subjecting 
the equations to a process referred to as “preconditioning.” This latter process is 
a device for replacing the system of equations by one mathematically equivalent 
but in which rounding off errors have a small effect on the final solution. 

The purpose of this paper is to describe in a very elementary way with simple 
examples the workings of the above mentioned phenomena. Some theorems 
which set out sharply the basis of the trouble will be stated and proved. 

General remarks. We begin with an illustrative numerical example: Consider 
the set of equations 


e+ty=1 
x+1.01y = 2. 


(1) 


These have as a solution y =100, x = —99. The solution is very sensitive to the 
value of the coefficient of y in the second equation. For example, if 1.01 is re- 
placed by 0.99, a change of 2%, the solution becomes y= —100, x=101, a 
tremendous change in the solution. 

We can replace the set of equations (1) by a second set which is mathe- 
matically equivalent but which is quite insensitive to changes in the coefficients. 
For example, if one subtracts 97 times the first equation of (1) from 100 times 
the second and then subtracts 700 times the second from 704 times the first one 
obtains the equations, 


3x + 4y = 103 


(2) 
4x + 3y = — 696. 


The set (2) is mathematically equivalent to the set (1) but is obviously less 
sensitive to small percent changes in the coefficients than in the set (1). One can 
say that replacing (1) by (2) is “preconditioning” the equations for computation. 
It must be clearly understood that preconditioning is only useful where the errors 
introduced by changes in the coefficients are at the control of the computer. For 
example, if the equations (1) are obtained from a physical experiment in which 
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all the coefficients are measured quantities and the coefficient 1.01, for instance, 
means a quantity between .99 and 1.03 then there is an essential indetermination 
in the numbers x and y which no amount of preconditioning of the equations 
can help. In this case, perhaps, all our abuse should be heaped upon the physical 
experiment rather than upon the equations. Where preconditioning is useful is 
in those cases where errors are introduced by the computing process itself. The 
most obvious example is the case where at each stage in the computation the 
numbers obtained are rounded off. 

The phenomenon of “ill condition” is more usually ascribed to a matrix than 
to a set of linear equations, the matrix in the case of a set of linear equations 
being that of the coefficients of the variables. In [1] and [2] of our bibliography, 
various measures of “ill condition” are given. The one due to von Neumann and 
Goldstine, namely the ratio of moduli of the largest to the smallest of the char- 
acteristic roots, will be given an interesting justification in Theorem 3 of this 
paper. In passing, we note that the von Neumann and Goldstine condition 
number of the matrix of equations (1) in our example is about 400, while that 
of equations (2) is 1. 

In general a system of linear equations with matrix A has a solution whose 
sensitivity to changes in the coefficients is roughly proportional to the elements 
of A-'. More explicitly, if the equations are 


= (¢ = 1, 2,3,--+,m) 
j=l 
and A =(a,;) with A~!=(c,;), then 


Ox; Ox; 
0b; O00 


Roughly speaking, an ill-conditioned matrix has for its inverse a matrix whose 
elements are relatively large. In the work that follows we will show that the 
computation of an approximate inverse to an ill-conditioned matrix gives rise 
to a complication, namely, that it is necessary to distinguish between an ap- 
proximate left inverse and an approximate right inverse. An example will be given 
to show that if linear equations are solved by matrix inversion, ill-condition is 
not as serious as might first appear provided one is careful not to confound 
approximate left inversion with approximate right inversion. 

Theoretical development. We first prove a theorem to show that there are 
matrices whose determinants are arbitrarily sensitive to changes in one of the 
elements: 


THEOREM 1. Let € and K be two arbitrary positive real numbers. For each n=2, 
there exists a non-singular (n Xn) matrix with the property that multiplication 
of a specified one of its non-zero elements by 1+ will multiply its determinant by 
1+K. 


c 


| 
| 
| 
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Proof. Let the matrix be A =(a;;) (¢, 7=1, 2, - - - , m); let the specified ele- 
ment be ai and let A;; be the cofactor of a;;. We are to choose A so that the 
following equation amongst determinants holds: 


+ €) Gin | Gin 
|= (1+K)] - 

Ani * Ann GQni Gn2*** Onn 


This is equivalent to where U=a2An+azAut - 
+4a:nAm. Here neither Ay, nor U involves the element ay. The value of an is 
given by 


—KU 


(3) 

— e) 
For K #e, we need only choose the elements of A (except for a) in such a way 
that U and Au are distinct from 0 and then choose ay by equation (3). If 
K =e the theorem is true for any matrix A of the form 


a1 0 0 ---0 

G22 Gon 
A=|- 

GQni * Onn 


provided A has a non-zero determinant. 


THEOREM 2. Let € and K be two arbitrary positive real numbers. For eachn2=2 
there exist non-singular (n Xn) matrices A and X such that XA has each of its 
elements differing from the corresponding element of the identity matrix by less 
than € and such that AX has all of its elements greater in absolute value than K. 


Proof. We will prove the 2 X2 case here leaving the general proof to the end 
of the paper. This is done because our next theorem concerns only 2X2 matrices 
and is a companion theorem to this one. 

Let e* be any number such that 0<e*<min (1, €) and let K* be a number 
such that K*>max (K, 2). 

Let 


| 

! 

| 

| 
— 
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and let 
K* + K* 
K*2 K*2 K* 
K* + e* +1 € K 


By direct computation, 


XA = ), 
1 — 


so that XA has its elements differing from those of the identity matrix by less 
than e. 
Also by direct computation 


(K** 2K* K*? ) 
+1—eé* +e ) 


é* 
K*2 2 K*2 2K*2 


Our choice of K* and e* assures that the magnitude of every element of AX is 
greater than K. 

The phenomenon exhibited in Theorem 2 is always associated with ill-condi- 
tioned matrices. In the next theorem we will show that the von Neumann and 
Goldstine condition number for such matrices is large (t.e. for small ¢ and large 
K). We will normalize the problem so that the characteristic roots have a 
product equal to unity in which case we may replace the ratio of the character- 
istic roots by the difference. 


THEOREM 3. Let ¢ and K be arbitrary positive numbers such that K/e>4. Let 
A and X be two (2X2) matrices with the properties: 

(1) A has non-negative elements and det A =1. 

(2) XA has all its elements differing from those of the identity matrix by e. 

(3) AX has all its elements greater in magnitude than K. 
If x; and x2 are the characteristic roots of A, then 


Proof. Let 


a b 


6 

\ 

| 

| 

| 
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a, b, c, d are all non-negative and ad—bc =1. We are given that 
tte te 
XA = ( 
ite 
The proof is similar for all changes of signs so we will consider only the case 


where 
1 
KA = T= ( 


€ 1—e 


Then X =TA~—!, so that 


—e(a + b)? ). 
e(c + d)? 1 — + b)(c + d) 
By hypothesis; 
ea+b)?> kK 
e(c+d*>K 


i+eda+d)(c+d>K 
ea+ b)(c+d)—1> K. 


From the first two of these inequalities 


K K 
o+b>4/— ct a> 4/— 
€ € 
Using ad—bc =1, the characteristic roots of A satisfy x*?—(a+d)x+1=0. This 


has roots x; and x: such that |x:—x2| =./(a+d)?—4. We also have c>\/K/e—d 
and K/e—a. 


If both \/K/e—a, \/K/e—d are positive we have 


(a + d) — + od bc 
€ € 


or 


or 


= 


| 
} 
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Hence 
€ K 
or 
K € 
(a+ >—+— +2. 
€ K 
Hence 
K € 
(a+d)*-4>—+4+—-2, 
€ K 
and 


K K 
€ K € 


If at least one of \/K/e—d, \/K/e—a is negative, suppose »/ K/e—a is negative. 
Then a?>(K/e). Hence, (a+d)?—4>a?—4>(K/e) —4 and, | x1 >VK/e—4. 

The author has proved a similar theorem under weakened hypotheses e.g., 
“XA has all of its elements differing from those of the identity matrix by 
e,” may be replaced by “XA has one of its elements differing from the cor- 
responding element of the identity matrix by e while the remaining elements of 
XA differ in absolute value from the corresponding elements of the identity 
matrix by an amount less than or equal to e.”’ Also the statement ‘‘AX has all its 
elements greater in magnitude than K” may be replaced by “AX has at least 
one of its elements greater in magnitude than K.” The proofs are considerably 
more complicated and would not add anything to the point of our discussion. 

Before giving a numerical application we may note the following consequence 
of Theorem 2. For a specified degree of accuracy an ill-conditioned matrix A 
may be associated with matrices X and Y where XA and AY differ from I by 
at most the specified error while AX and YA may differ considerably from J. 
We may then say that X is a suitable left inverse but not a suitable right inverse 
of A while Y is a suitable right inverse but an unsuitable left inverse of A. If 
one is to solve a system of linear equations by matrix inversion it seems reason- 
able that one should not confuse X and Y. More specifically, consider the set of 
equations 


= b; (i = 1,2,3-+-n). 


j=l 


As a matrix equation this may be written 


AU=B8B where A = and 


| 

| 

| 
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uy b; 
B= bs 
bn 


Here a solution is U=A~'B. If we are approximating to A~! it is obvious that 
what is required is an approximate /eft inverse of A. On the other hand if we 
consider the equations 


= 6; (j = 1,2,+++ m) 
t=—1 
and write them in matrix form VA=C, where A =(a;;), V=(01, v2, - + * 5 Ua), 
C=(c1, C2, Cn), the solution V=CA™-! requires that A~! be approximated 


by an approximate right inverse of A. A numerical example which follows il- 
lustrates the error introduced by confusing an approximate left inverse with an 
approximate right inverse. 

Numerical example. Consider the set of equations 


x+y+2= 6.00 
(4) y+ 1.012 = 6.03 
x + 0.99y + 2 = 5.98. 


These equations have an exact solution x =1, y=2, z=3, Equations (4) can be 
written in the matrix form: 


1 1 1 x 6.00 
(5) 1 1 1.01|| y| = |6.03). 
10.99 1 z 5.98 
The matrix 
1 
A=/1i1 1 1.01 
1 0.99 1 


is easily seen to be an ill-conditioned matrix. If we do not invert exactly it ap- 
pears that an approximate left inverse is more appropriate than an approximate 
right inverse. We will assume computations are carried out in some systematic 
manner and rounded off to two decimal places. Conceivably, one might com- 
pute as an approximate left inverse the matrix 
1.01 —100 100 
X= 99.99 2 102}. 


— 162.01 102 0 


= 
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Here 
1.01 .01 
XA ={|-.01 1.01 .01 
—.01 —.01 
while 
—2.01 4 —2 
AX = |-2.03 5.02 
—2.01 3.98 —0.98 


Hence X is a good approximation as a left inverse of A but a poor approximation 
as a right inverse. If we use X to solve (5) we obtain 


6.00 
y| = X|6.03 
5.98 
which yields x =1.06, y=2.04, z=3.00, which is a good approximation to the 


exact solution considering the condition of A. 
In the same way the matrix 


1.01 —101.99 100.01 


Y= 102 2 — 100 
102 100 0 
has the properties 
1.01 
AY =|-—.01 1.01 .01 
—.01 —.01 1.01 


and 
—0.97 -—1.97 —1.99 
YA = + 5 4.02 
—2 —2 —1 
so that Y is a good approximate right inverse but a poor approximate left inverse. 


If one ignores the distinction between an approximate left inverse and an 
approximate right inverse and tries to use Y to solve (5) in the form 


6.00 
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one obtains x = — 10.88, y =26.06, z= —9 an answer with a fantastically large 
error. 

Since, for computing purposes, distinction between left and right approximate 
inverse is of some importance (and of utmost importance in the case of very 
badly conditioned matrices) it is worthwhile to look into some of the methods 
of matrix inversion to determine what kind of inverse they determine when used 
approximately. For example; Albert in [5], gives a method which when used 
approximately gives a left inverse while Andree in [6] gives a method which 
approximates to a “hybrid” inverse. It is not easy to see in most cases whether 
the method of approximation leads to a one-sided or hybrid inverse. Other 
things being equal a one-sided inverse is preferable to a hybrid. Where iterative 
processes are used it is advisable to calculate as a left inverse and check each 
stage for accuracy as a right inverse. In the case of symmetric matrices any 
method of inverse approximation which yields a symmetric matrix as an inverse 
will be good since the resultant matrix has the same error both as a left and as a 
right inverse. These remarks are mainly of use because one cannot tell in many 
cases whether or not a matrix is ill-conditioned. Many innocent looking matrices 
are very badly conditioned. Finally, where a set of equations has an obviously 
ill conditioned matrix and preconditioning can easily be carried out, it is ad- 
visable to do so. 

We close with a proof of Theorem 2 in the general case. 

Proof of Theorem 2 for (Xm) matrices. Let e* be any number such that 


0<e*<min (1, €). Let A =(a,;) (¢, 7=1, 2, - +--+, m), where a,; are to be deter- 
mined later in such a way that det A =1. 
Let A-'=A,; (i, 7=1, 2, --+,m). Define X by the equation 
1+ ¢, 
1+ é, At. 


Then X =(J+£)A-' where J is the identity matrix and 


* * 

= 

* * 

&, #&, ,€ 


Denote E by (li;) where (i, j=1, 2,---, m). It follows that AX 
=I+AEA—!=I+Y where Y=AEA~—. If Y=(y,,;) (¢,7=1, 2, +>, we have 


= > GivbuvA = ain) ( Au). 


tan] ve] 


| 
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Suppose we can choose A so that 

(1) det A =1. 

(2) The sum of the elements in each row of A is greater than 2K/e*. 

(3) All sums >>", A,; are non-zero integers. 
Then all the y,; are greater in magnitude than 2K, so that all the elements of 
AX are greater in magnitude than 2K —1. Let K* be a positive integer satisfying 
the conditions K*>2 and K*>2K/e* and put 


1 K*---K 


Baul. 
0 0 O---1 

Let the inverse of B have elements B,; (¢, 7=1, 2, 3, -- +, m). Then B has the 


properties (1) and (2). With regard to property (3) we can say that all sums 
>"*., B,, are integers at least one of which is different from 0. For if all 
>-*., Br, were equal to 0, suppose we choose Z so that ZB=I+E. Then BZ 
=J+BEB-', We have BEB-!=0 by the condition >>*., B,,=0 for all s. Hence 
Z is an exact right inverse of B, a contradiction, since ZB =I+E and E¥0. The 
matrices B and Z then have the properties: 

(a) Each element of ZB differs from the corresponding element of the 

identity matrix by an amount less than e. 
(b) Not all elements which are off the main diagonal of BZ are equal to 0, 
and each such non-zero element of BZ is greater in magnitude than K. 

If BZ has no element equal to 0, it is also true that every element, including 
those on the main diagonal, is greater in magnitude than K. In this case B 
satisfies the conditions imposed on A by the theorem. 

If BZ has at least one zero we proceed as follows: Let Cs, C3, ---,C, be 
indeterminates to be replaced by positive integers. Form a matrix A as follows: 
the first row of A is taken to be the same as the first row of B. For 
(¢=2,3,---,m), ({=1,2,3,+--,m), take Then A satisfies the 
required conditions (1) and (2). With regard to the condition (3), each sum 
>-*., Are is a linear combination of the By with coefficients products of the 
C; and such that at least one of the B,; in the sum is different from 0. Hence, 
for each s, ; ae A,, is a polynomial in C2, C3, - - - , C, which is not identically 
0. As such a polynomial equated to 0 has at most a finite number of solutions 
in the C;, there exist choices of C2, C3, - - - , C, such that each C; is a positive 
integer and such that for each s, p> A,,~0. For such a choice the matrix A 
satisfies the conditions of the theorem. 
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CONSISTENCY AND COMPLETENESS—A RESUME 
FRANK DeSUA, Ohio University 


The modern foundational investigations of mathematics, encroaching as they 
do on unresolved philosophical questions, have had their polemic aspects. Con- 
troversy over ontological problems gave rise at the turn of the century to the 
three “schools” associated with the programs of Russell and Whitehead, of 
Hilbert, and of Brouwer, a trichotomy which is preserved to some extent even 
today. However, those results of the investigations which are perhaps the most 
significant for the mathematician, namely those embodied in the famous in- 
completeness theorems of Gédel, must be considered as well established and 
unimpeachable as any found in the whole of mathematics, based as they are on 
methods acceptable to the devotees of each of these programs. 

Among the values to be gained from these investigations is the insight 
afforded into the nature of our mathematical preconceptions, and the revelation 
that even some of the most compelling of these may actually be misconceptions. 
Of course, once this is discovered, the void created by the dissolution of these 
misconceptions is apt to be immediately filled by newly formed assumptions, 
consciously or otherwise. The history of mathematics, as of any science, is to 
some extent the story of the continual replacement of one set of misconceptions 
by another. This is of course no cause for despair, for the newly instated as- 
sumptions very often possess the merit of being closer approximations to truth 
than those which they replace. 

We are concerned in this paper with the detection of one such misconception 
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—a misconception concerning the scope and power of formal mathematical 
methods—and Gédel’s remarkable discovery of largely unsuspected limitations 
inherent in such methods. This discovery is part of, and a product of, the “crisis 
in the foundations” of which the rumblings have not yet died away. It is not 
yet entirely clear what new assumptions are to supplant the old, but it is cer- 
tainly clear that the old must give way. 

Our paper is divided into three sections, the first of which is entitled “the 
Problem,” which is the problem of the consistency and completeness of mathe- 
matics, first explicitly formulated by Hilbert. The last section “the Solution,” 
concerns the surprisingly negative solution of Gédel. The second section, “For- 
mal Systems,” is a digression necessary for an appreciation of the nature of the 
solution. Before we attempt a statement of the problem, we concern ourselves 
to some extent with how it arose. 

The problem. The theory of abstract sets or classes may be said to be the 
basic discipline of pure mathematics. Belief in the existence of sets is almost 
universal among mathematicians, a belief which seems to be preordained by the 
dictates of common sense. So compelling is the notion of set that even those 
who harbor doubts must at.least feign belief part of the time in order to dis- 
course intelligibly. Yet this concept, despite its seeming necessity and superficial 
innocence, has proved to be the Pandora’s Box of mathematics. 

Consider how sets are specified. Two methods are conceivable. A set may be 
specified extensionally, i.e., by a listing of the members such as “the set consist- 
ing of the numbers 2, 3, 5, and 7” or intensionally through some defining prop- 
erty or necessary and sufficient condition for membership—technically known 
as a “monadic propositional function”—as “the set of all primes less than 10.” 
The extensional specification of sets is, of course, precluded as soon as we con- 
sider infinite sets. 

The intensional method is often canonized in the form of a so-called logical 
principle, the (unrestricted) Principle of Abstraction, which states that every 
property determines a set, or in technical language, every monadic propositional 
function has exactly one set as extension. Symbolically, this becomes: 


[2 P(x)] 


where S is the set determined and P is the defining property. 

Briefly the argument in support of this principle is this: Every condition or 
property determines the set of all entities possessing the property or obeying 
the condition (even in the case where no such entity exists, the set determined 
being empty). Conversely every set determines at least one property, namely 
the property of being a member of the set, if no other. 

The principle is widely appealed to in mathematics—often in some restricted 
form—either explicitly or tacitly. It is implicit, for example, in the high school 
geometry definition of a “locus” as a point set every element of which obeys 
some metric condition. 

In the principle of abstraction, then, we have an intuitively obvious prin- 
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ciple of wide applicability. Unfortunately there is one thing wrong with it. In 
the unrestricted form stated above, it happens to be false. It was proved to be 
false by Bertrand Russell in 1901. 

Let’s examine the process of set formation as embodied in the principle of 
abstraction. Out of the aggregate of all entities we choose those which possess a 
given property to form our set. The flaw in the principle of abstraction lies in 
the fact that it ignores one entity, namely the newly formed set. Conceivably 
the property used to form the set may neither apply nor fail to apply to this set. 
In this case no set is actually determined. Russell discovered a property, pos- 
sessed by many entities, yet which cannot consistently be said to determine any 
set. 

The property in question is the property which a set x may possess of “not 
being a member of itself” (x@x). Many sets possess the property—a set of 
material objects, for example—while other sets fail to possess the property—the 
set of all sets would, assuming the existence of such a set. Indeed in an unre- 
stricted universe of discourse (assuming the allowability of such a universe), 
there must be as many sets of one kind as of the other, for the complement of 
any set possessing the property is a set which does not, and conversely. 

It is easily seen that the above mentioned property determines no set, for 
suppose there were such a set S. From the defining property, xGSex€x. 
Does S itself possess the property, or not? Putting S for x we have, SES 
«SES, which is a contradiction. Therefore no such set S exists, although 
innumerable sets x possess the property x€x. 

Actually, the first proof of the non-validity of the principle of abstraction 
came several years earlier in a result by Burali-Forti concerning ordinal num- 
bers. Sets of numbers of various kinds are common objects of study in mathe- 
matics. It would seem then that the set of all ordinals is perfectly well defined 
according to the canons of common sense mathematics. Given any entity we 
can decide whether it is a member of the set or not if we can determine whether 
it is an ordinal number or not. It came therefore as a surprise when in 1897 
Burali-Forti demonstrated the non-existence of the set of all ordinals. 

To recapitulate the proof, recall the following theorems of ordinal number 
theory: Any initial segment of the ordinals, that is, any set of ordinals contain- 
ing all ordinals preceding any of its members, can be ordered according to the 
magnitude of its members. The resulting ordered set is well-ordered; conse- 
quently it has an ordinal number. Further this number exceeds any member of 
the set (in fact, it is the least number which exceeds every member of the set). 

Now suppose the property of being an ordinal number determined a set S, 
t.e., let S be the set of all ordinals. Then S would be a set to which the above 
theorems apply, so that it would have an ordinal number Q. Accordingly, Q 
would exceed every member of S. But if S is the set of all ordinals, 2 must be 
one of its members. Consequently, 2 exceeds itself, which is a contradiction. 

Superficially somewhat similar difficulties arise, not only in the theory of 
sets, but also in the theory of statements or propositions. It would seem that 
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any statement which purports to be a statement of fact must be either true or 
false. If the statement conforms with the facts it is true, otherwise false. Thus, 
for example, the statement “John has quit smoking” is true only if John has 
actually quit smoking. (There are difficulties of an existential nature in some 
cases, as for example with such statements as “Pegasus has wings” or “Four 
sided triangles are isosceles.” But such difficulties are capable of resolution 
through the adoption of certain conventions of interpretation, as in Russell’s 
and Quine’s theories of singular terms and descriptions.) 

This view of the validity of statements, however, is subject to difficulties of 
another kind if no restrictions are imposed on statement formation. For if we 
are talking about statements themselves, 7.e., if we consider statements about 
statements, as we often do in mathematics, the statements themselves become 
part of the facts. Conceivably in an unrestricted discussion the facts may be so 
affected as to render certain such statements about statements neither true nor 
false. 

That this is actually the case was discovered by the ancient Greeks. They 
constructed a statement (the Epimenides), of the form “This statement is 
false,” which asserts its own. falsehood. This statement, referring as it does to 
itself, can only be true if it is indeed false, and conversely. 

It might be pointed out that such Epimenidean statements cannot be con- 
structed in any properly restricted mathematical system, so that the contradic- 
tion does not arise in practice. (The Russell and Burali-Forti contradictions 
actually did arise in mathematical discourse, and still do occasionally.) How- 
ever, as we shall see, something very closely resembling such Epimenidean state- 
ments actually can be constructed in certain formal systems, namely Gédelian 
formulas. This is a profound syntactical fact which we shall examine later. 

A great many of these demonstrations running counter to intuitive prin- 
ciples have unexpectedly come to light. Cantor proved the non-existence of the 
set of all cardinal numbers. Berry and Richard have shown the non-existence 
of certain apparently well defined denumerable sets of finite numbers. Skolem 
has shown conversely the existence of certain number sets possessing no defining 
property. 

We have taken the position here that these “contradictions” actually con- 
stitute demonstrations, by reductio ad absurdum, of the non-validity of certain 
intuitive principles of long standing, such as the principle of abstraction. Other 
positions concerning the significance of these contradictions are widely held. 
The authors of the contradictions themselves very often tacitly defended these 
principles, and consequently regarded the contradictions as “paradoxes.” They 
had then to revise their concept of “property” or “condition,” rather than “set,” 
etc. Whatever position is taken however, it is clear that unbridled intuitive rea- 
soning, even on the most basic level, can and sometimes does lead to difficulties. 

With this realization we are now in a position to state what was until re- 
cently a central problem of the modern foundational investigations of mathe- 
matics, namely to devise restrictions on intuitive reasoning sufficiently strong 
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to preclude the possibility of inconsistency but sufficiently weak to aliow mathe- 
matics the greatest possible latitude for development, at least to allow the 
development it has already undergone. Further the problem is to prove, if 
possible, by unimpeachable methods (and therein lies a danger—the principles 
discredited by Russell were “unimpeachable” in their day) that these restric- 
tions are actually sufficiently strong, i.e., that the resulting mathematical sys- 
tems are consistent. 

With the first part of this problem we shall not concern ourselves except to 
say that restrictions have been devised which seem to fulfill the requirements. 
There are Russell’s type theory, Ramsey’s simplification of this theory, Quine’s 
stratification, and Zermelo’s axiomatic set theory to mention a few. The ade- 
quacy of the resulting system for the derivation of classical mathematics is 
proved very directly. Either classical mathematics is actually derived within 
the system (as in Russell and Whitehead’s Principia Mathematica) or it is shown 
that any theorem that can be derived in, say, Principia, can also be derived in 
the system under consideration. The proof of the consistency of these systems 
is quite another matter, and this is the part of the problem with which we are 
mainly concerned. 

Formal systems. The systems whose consistency is investigated are not 
actual systems of ordinary intuitive mathematics but rather certain specializa- 
tions thereof known as formal systems or formal languages. We shall investigate 
this concept in some detail presently, but here we want to point out differences 
between formal mathematics and the postulate systems of ordinary mathe- 
matics. In ordinary mathematics it is helpful in carrying out deductions to know 
the meanings of the technical terms (or “undefined” terms) involved. The se- 
mantics of the system is of consequence. Even when the system is semi-formalized 
with regard to technical terms, so that they become truly undefined, semantics 
plays a heavy role in the logic involved in the deductions, which is usually on a 
primitive verbal level. 

In formal mathematics on the other hand the technical terms are analyzed 
to the point where all their properties which matter for the deduction of theorems 
have been expressed by axioms. The terms must then have a recognizable form 
rather than an apprehensible meaning. It is then possible to treat these terms 
as words essentially without meaning, for to say that these terms have meanings 
essential to the deductions amounts to saying that not all their properties have 
been expressed by axioms, i.e., that the analysis is incomplete. Even the logic 
of the system is explicitly incorporated into the system in the form of axioms 
and rules of inference, and the whole is rendered in symbolic rather than verbal 
form. When this analysis is completed, and meanings thus left out of account, 
questions of form or structure become the sole consideration. Semantics has 
been replaced by syntax. 

There are several good reasons for this preference for formal over intuitive 
systems in many foundational investigations, not the least of which is the fact 
that intuitive mathematics is a function of the mathematician as well as of the 
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subject matter. Whereas a topologist may embrace such principles as the axiom 
of choice or the continuum hypothesis without hesitancy, an algebraist may do 
so with compunction. Thus it becomes of interest to compare systems in which 
these principles are explicitly included with systems in which they are excluded. 
Also standards of rigor vary with the times and the mathematician, one man’s 
proof being another man’s sophistry. The intuitive mathematics of the great 
eighteenth century analysts, for example, is demonstrably inconsistent. Ques- 
tions of rigor do not arise however in a system restricted to finitary and con- 
structive methods of proof as is the case in many formal systems. Perhaps the 
chief reason for the use of formal systems rather than intuitive systems, how- 
ever, is that only within formal systems can such vague semantical concepts as 
axiom, theorem, proof, consistency, etc., be made explicit and precise. Formaliza- 
tion is a necessary concomitant to precision. 

The study of formal systems has been criticized (and extolled!) on the 
grounds that it “reduces mathematics to a meaningless game.” It seems plausible 
that such criticism is based on a misconception concerning the purpose of these 
studies. It would be as much to the point to criticize the study of topological 
spaces on the grounds that such “spaces” are incapable of propagating electro- 
magnetic radiation. It must be kept in mind that the partial formalization of 
much of mathematics—of algebra, for example—has proceeded steadily since 
the time of the Babylonians, yet algebra certainly remains meaningful mathe- 
matics. Abstraction of meaning does not imply loss of meaning. 

The extent to which intuitive mathematics can be formalized is a problem 
on which much light has been shed by recent foundational investigations. That 
at least a very considerable measure of formalization is possible for mathemati- 
cal theories is evident upon a perusal of some of the attempts at such formaliza- 
tion—Whitehead and Russell’s Principia Mathematica, for example. But that 
limitations do exist on the kinds of problems which can be stated and proved 
within a formal system is revealed by the researches of Gédel, Church, and 
others. And the delineation of such limitations is the subject of continuing 
research. Involved in the process of formalizing a theory may be difficulties as 
great as those involved in the subject matter of the theory itself. 

Let’s examine the way in which a formal system is actually set up. We start - 
with a certain set of symbols which may be infinite in number. Any finite linear 
array of these, allowing repetitions, is called a formula. The formulas are di- 
vided into two disjoint sets called well-formed formulas (w.f.f.) and gibberish 
(only the w.f.f. are to have meaning in a semantic interpretation of the system). 
The w.f.f. are defined recursively, and in a purely syntactical or morphological 
manner, 7.e., solely in terms of their structure. Moreover they are defined effec- 
tively, t.e., a finite process is specified whereby we can determine whether a 
given formula is well-formed or not. 

A certain subset of the w.f.f. are chosen to be our initial axioms. (The age- 
old criterion of “self-evidence” of axioms is dispensed with, there being no de- 
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grees of perspicuity in a system from which all meaning has been abstracted.) 
New w.f.f. are obtained or “derived” from these axioms by the application of 
certain rules of inference, which again are recursive, syntactical, and effective. 
By a proof we mean a finite sequence of w.f.f.’s, each of which is either an axiom 
or is obtainable from one or more earlier formulas of the sequence by an ap- 
plication of a rule of inference. 

The last formula of a proof is, by definition, a theorem. We write A to 
denote that A is a theorem. (The symbol ‘t+’ is called the assertion symbol.) 
Whether a given sequence of formulas constitutes a proof or not can be effectively 
determined, but whether a given formula is a theorem or not generally can not 
be determined by any effective uniform procedure. 

An actual example of a formal system may serve to clarify some of these 
definitions. The following is one possible formalization of a fragment of modern 
logic known as the “propositional calculus”: 


Symbols: An infinite list of lower case letters, primed and unprimed: a, },c,---, 
a’, b’, c’, a”, b”, c”, 
parentheses: (, and ) 
negation symbol (tilde): ~ 
alternation symbol: V 
Any finite array of these, for example \/))a’~, is a formula. This 
particular one happens to be gibberish. 
wff.: Any letter (primed or unprimed) alone is w.f. If A and B are w.f., so 
are ~A and (A\/B). Examples: r is w.f. So are ~7r and (r\V~7). 
Axioms: An infinite list of axioms which come under the following axiom 
schemata: 
(1) (AV~(AVA)) 
(2) ((~BVC)VB) 
(3) (~(DV ~E)V(~(EVF)V(FVD))) 
Each of these schemata represents an infinite set of axioms, namely one 
for each possible choice of w.f.f. in place of A, B, C, D, E, and F. 
Rule of inference: From two formulas of the form G and (~G\/), we infer the 
formula H. 

This system is capable of a semantical interpretation in terms of the so- 
called logical connectives “not,” “either - - - or - - - ,” etc. The letters may stand 
for sentences or propositions. Thus if we let “r” denote “it is raining,” then 
“~r” becomes “it is not raining” and “(r\/~r)” becomes “either it is raining 
or it is not.” As an example of a theorem and its proof we assert that F (r\/~7). 
Proof: 

(i) 
[from axiom schemata 1 replacing A by ~r] 

[from 3 replacing D by ~r, E by (~r\/~1), and F by r] 
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[using the rule of inference on steps (i) and (ii) ] 

(iv) ((~rV~r) Vr) 

[from 2 replacing B by r, and C by ~r] 
(v) 
[using the rule of inference on steps (iii) and (iv) ] 

Thus we have “proved” that it is either raining or it is not. 

We note here, for later use, a binary connective of considerable importance, 
namely the conditional or material implication, symbolized by “—,” and trans- 
lated by one sense of the idiom “if --- then ---.” If A and B are w-f., then 
(A—B) is w.f. A rule of inference equivalent to that given in the above system, 
but much more widely used, is modus ponens (or strictly speaking, its syntacti- 
cal analogue, detachment) which allows us to obtain from two formulas of the 
form A and (A—B) the formula B. 

Two further definitions are important for our discussion. A formal system is 
consistent if there is no w.f.f. A such that both FA and /~A. (In what follows 
a stronger condition called w-consistency is sometimes required, but we shall 
forego an explanation of this for the sake of brevity.) A formal system is com- 
plete if, for every w.f.f. A, either A or the system would become inconsistent 
if A were introduced as an additional formal axiom (the reader may think of 
this as either FA or F ~A, 1.e., the system is complete whenever every w.f.f. 
is theoretically “decidable,” although strictly speaking the two concepts do not 
coincide in all formal systems). 

The solution. Our problem now becomes the problem of setting up a formal 
system which is adequate for the derivation of mathematics, and which further 
is consistent (1.e., leads to no contradiction) and complete (i.e., we would like 
every mathematical proposition to be at least theoretically decidable in the 
sense that it can either be proved or refuted). The problem of adequacy was 
solved by Frege, Whitehead, and Russell. (Of course, the adequacy of their 
systems was demonstrated for extant mathematics only. Nothing can be ascer- 
tained concerning the adequacy for all future developments of mathematics.) 
Several systems have been proposed which it is hoped are consistent as well as 
adequate. To mention a few, there is Russell and Whitehead’s Principia Mathe- 
matica, Quine’s Mathematical Logic, and systems by Bernays, Fraenkel, Gédel, 
von Neumann, Zermelo and others. Although it is hoped that these systems are 
consistent, none have actually been proved to be so. 

The consistency problem was partially reduced in classical mathematics by 
the so-called “method of models.” Thus Beltrami in 1868 reduced the con- 
sistency of Lobatchevsky’s non-Euclidean geometry to that of Euclidean geom- 
etry. In this and other manners, the consistency of large parts of classical mathe- 
matics is ultimately reducible to that of the arithmetic of natural numbers, 
particularly as this theory is embodied in the Peano axioms, or in a set theory 
sufficiently comprehensive for the derivation of the Peano axioms. 

For several inadequate partial systems the consistency problem has been 
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successfully solved. Post in 1921 demonstrated both the consistency and the 
completeness (in a somewhat different sense than that defined above) of the 
propositional calculus. The first order functional calculus was proved consistent 
in 1928 by Hilbert and Ackermann, and complete (in still another sense) in 
1930 by Gédel. More recently Tarski has demonstrated the consistency and 
completeness of the arithmetic of real numbers, and also of elementary geom- 
etry. 

With this background of success it was widely felt that the solution of the 
consistency problem for the arithmetic of natural numbers was essentially at 
hand, that it was now merely a matter of detail. The necessary tools had already 
been forged in Hilbert’s new branch of research, metamathematics, and few 
doubted that a positive solution to the problem was already within reach. Sur- 
prisingly however, despite intensive efforts by Hilbert and others, and very close 
approaches to a solution, the problem remained refractory. 

This puzzling situation was illuminated in 1931 by the appearance of certain 
remarkable and unexpected results of Gédel which solved in a negative fashion 
both the consistency and the completeness problem. These theorems, among the 
most remarkable in all of mathematics, revealed an unforeseen limitation in the 
methods of formal mathematics in some ways analogous to the limitations on 
physical measurements imposed by the Heisenberg Uncertainty Principle. 
Working within systems either presupposing the Peano axioms or possessing a 
syntax sufficiently rich for the derivation of these axioms—systems related to 
the Principia Mathematica—Gédel showed the impossibility of proving the 
consistency of such systems by so-called “finitary” methods, that is by construc- 
tive methods formalizable within the system itself. So the proof of consistency 
which seemed so close to attainment is actually impossible of attainment, at 
least in the form originally envisioned by Hilbert. 

This remarkable result is but a consequence of even more fundamental 
results concerning the completeness of the formal systems studied by Gédel. 
Gédel demonstrated the incompleteness of such systems, that is, he revealed the 
existence of “undecidable” propositions, of which consistency is one. This result 
is particularly surprising in view of the fact that the formal systems with which 
Gédel worked are sufficiently rich in syntax for the derivation of all of classical 
mathematics (and presumably of much of mathematics yet to be developed). 
Yet despite this richness of syntax, it is possible to formulate within these sys- 
tems arithmetical problems which are theoretically unsolvable within the system. 

No attempt will be made here to outline Gédel’s demonstrations, for even 
such an outline would of necessity be long, detailed, and technical. Instead we 
will content ourselves with a heuristic account of Gédel’s reasoning which we 
hope will at least reveal the “point” of the proof. 

Suppose we denote by “CONSIS” a syntactical formula whose semantic 
interpretation implies the consistency of the formal system (for example, 
CONSIS may mean “there is no proof for 1=0”). The consistency problem 
would then be related to the problem of demonstrating / CONSIS. One pos- 
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sible line of attack on the problem would be to seek a sufficient condition A 
such that (A—*>CONSIS) and then to try to demonstrate -}A. CONSIS 
would then follow by modus ponens. 

Now Gédel foresaw what others had failed to suspect, namely that CONSIS 
was neither provable nor refutable if the system is truly consistent. To demon- 
strate this he hit upon the following ingenious scheme. Suppose we seek a condi- 
tion B, necessary for CONSIS [so that + (CONSIS-—B) | and then show that 
B is true but unprovable within the system. It would then follow that CONSIS 
is undecidable, for if H CONSIS, then B would be provable by medus ponens, 
which would be a contradiction. 

The question now remains, can we construct such a formal statement, or 
formula B which, assuming consistency, is true but unprovable within the 
formal system? If we could, we will have solved in a negative fashion both the 
consistency problem and the completeness problem, 7.e., we will have shown 
that the methods available in a formal system adequate for extant mathematics 
are inadequate to demonstrate its own consistency, and thus inadequate to de- 
cide even theoretically all mathematical propositions. 

Taking a hint from the Epimenides antinomy, Gédel constructed such an 
undecidable formula B by constructing a formula whose semantical content is 
that B is undecidable. Consider the Epimenidean statement: “This statement is 
false.” By analogy, we might try “This statement can neither be proved nor dis- 
proved” for our statement B. However, one difficulty in following this Epi- 
menidean construction in the construction of our Gédelian statement B is that 
the self-reference occurring in the Epimenidean statement is semantical in char- 
acter. Such semantical self-reference is of course precluded in a formal system. 
We must seek a syntactical device which will enable statement B to refer to itself 
in a purely structural manner. 

A semi-syntactical device is available in the English language through the 
use of quotation marks. The following statement, for example, is Gédelian in 
character: 

“yields a statement which can neither be proved nor disproved when 
appended to its own quotation” yields a statement which can neither be 
proved nor disproved when appended to its own quotation. 

Of course the device of quotation is a part of the syntax of natural languages 
and not of mathematics, formal or otherwise, so it was necessary to seek further. 
Gédel accomplished his aim by a very elaborate and ingenious numerical device, 
assigning coordinates to the decidable formulas of the formal system, establish- 
ing an enumeration of these formulas, and constructing the undecidable formula 
B by an application of Cantor’s diagonal process. 

These results on incompleteness refer of course only to certain formal sys- 
tems whose syntax is adequate for the formalization of many properties of the 
natural numbers used in the construction of the “undecidable” proposition B 
(of the Chinese remainder theorem, and of the possibility of mapping integers 
onto finite sequences of integers, for example). It might be thought that the 
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defect of the incompleteness of such systems could be remedied by the addition 
of further effective lists of axioms which would enable one to derive the proposi- 
tions which were previously undecidable. But very surprisingly even this ex- 
pedient, for all its apparent reasonableness, fails utterly. For Gédel has further 
shown that not only are the formal systems incomplete but also they are in- 
completable in the sense that new undecidable propositions can always be con- 
structed within the system enlarged by the addition of these new axioms. 

A final word concerning the significance of these syntactical results for ordi- 
nary intuitive mathematics. Gédel has shown that the formal systems known to 
be adequate for the derivation of mathematics are unsafe in the sense that their 
consistency cannot be demonstrated by finitary methods formalizable within the 
system, whereas any system known to be safe in this sense is inadequate. We 
are left in the position of the student of elementary mathematics who “knows 
what a logarithm is, but just can’t put it into words.” We know, or at least 
strongly suspect, mathematics to be consistent. But we just cannot incorporate 
this knowledge into a formal system. 

The methods used in Gédel’s demonstrations were constructive and finitary 
in character, and consequently acceptable even to the most cautious of the 
critics of mathematical methods. In view of Gédel’s results, if a proof of con- 
sistency is to be forthcoming it must employ methods which are not formalizable 
within the system. Such proofs have been devised. Thus Gentzen has proved 
the consistency of number theory using transfinite induction up to a sufficiently 
large ordinal—a method of proof not formalizable within the system whose con- 
sistency it proves. While such methods are perhaps acceptable to most mathe- 
maticians, the acceptance is by no means universal. Indeed there are those to 
whom no proof will be acceptable except one of a type explicity ruled out by 
Gédel’s discoveries. In view of the lack of a universally acceptable proof, belief 
in the consistency of mathematics becomes then somewhat a matter of faith 
rather than reason. 

Suppose we loosely define a religion as any discipline whose foundations rest 
on an element of faith, irrespective of any element of reason which may be 
present. Quantum mechanics for example would be a religion under this defini- 
tion. But mathematics would hold the unique position of being the only branch 


of theology possessing a rigorous demonstration of the fact that it should be so 
classified. 
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A CHAIN OF CIRCLES ASSOCIATED WITH THE n-LINE 
J. W. CLAWSON, Ursinus College 


1. Introduction. The Clifford-de Longchamps-Morley-Bath series of theo- 
rems in the geometry of the -line are well known and were stated in an earlier 
paper.* We now add another chain of circles to the system. In section 2 an 
analytic method is set up and applied to the Morley circles. In sections 3 and 4 
the de Longchamps points and the Clifford points and circles are discussed. 
In sections 5 and 6 a chain of circles which we shall call the Miquel circles and 
a chain of points which we shall call the Hermes points are defined. In section 7 
the new chain of circles is reached and its chief properties proved. In the final 
sections the special cases of the 5-line and the 6-line are considered. 

2. The Morley circles. | We take m lines in the plane, no three concurrent 
and no two parallel, and consider the “normal” equation of /;: x cos ay+y sin ay 
=p;. If z=x+iy and 3=x—iy, this can be put into the form 2+f?2 =k, where 
f is an abbreviation for cos ay+# sin ay, and ky is twice py. Throughout this paper 
the boldface numerals 1, 2,*--, m are used in this sense. For example 
1*=cos 3a,+7 sin 3ay. 

Thus /, is 2+1°92=1 and is 2+27%2=k,2. Then Ay, the intersection of /; 
and /, is found to be z= —12{ki2—kel }/(12), where (12) is an abbreviation for 
1?—2?. We shall use (123) as an abbreviation for (12)(13)(23), (2346) for 
(23) (24) (26) (34)( 36)(46), and so on. The numerals will always be placed in 
ascending order. 

If, asin this ase, z is a function of k’s and boldface numerals, 2 can be found 
by replacing 1, 2,---, m by 1/1, 1/2,--+-+, 1/mn, since cos ay—i sin ay 
=1/(cos sin ay). Thus, for Ai, 2= {kil —h:2}/(12). 

The square of the distance between two points (x1, y1) and (x2, ye) is given by 
(21 — 22) (1 —Z2). Hence the equation of the perpendicular bisector of A12A23 is 


This reduces to 


~ (12)(23) 
Similarly the perpendicular bisector of A13A¢3 is 


2% 


[4112(23) — — 1232} — &,23(12)]. 


* This MonrTBLY, vol. 61, 1954, p. 161. 

t These circles were first announced by de Longchamps, Nouv. Corres. de Math., 1877, p. 306. 
Because of an important paper by F. Morley, Trans. Amer. Math. Soc., 1900, p. 99, they have been 
called Morley circles. R. Goormaghtigh in a paper in this MONTHLY, vol. 54, 1947, p. 327, calls them 
the centric circles. Other references are found in H. F, Baker, Introduction to Pure Geom., note 1, 
p. 330, Cambridge, 1943, 
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—3 
(13) (23) 


Eliminating 2, we find that Cis3;, the center of the circle circumscribing 


A2A13A23 is 


[2113(23) + + ks {3* — ]. 


23 yee. — k213(13) + 
(123) 
The radius of this circle is 
k 
123 { 11(23) — ke2(13) + 
(123) 


The equation of the circle is 
(123)2z — — + 
— 2123{ k,23(23) — k213(13) + h312(12)} 
= 123{ kike3(12) — kiks2(13) + keksl(23)}. 


We shall abbreviate these expressions so that Ciss is 123 }&:23(23)/(123), riss is 
123 >>k1(23)/(123), and the circle is 


(123)2z — z >> hil*(23) — 2123 >> = 123 by ks3(12). 


In these summations the numerals in the subscripts of the k’s and in the 
boldface numbers are to be rearranged in every possible way except that they 
must always be in ascending order among the k’s and the boldface numbers re- 
spectively; the sign of a term is plus or minus according as the number of trans- 
positions needed to put the numerals back into the original order is even or odd. 

We shall use n! for 12 - - - n, n!/f for n! with f omitted, n!/fg for n! with 
f and g omitted, n!/f! for f+-1 £+2 - - -n and so on. The same abbreviations 
will be used for integers. 

It may be noted that (n!/12)(n!) = (m!/1)(m!/2)(12). We now set up a chain 
of points as follows. Cy! is defined as (—1)*—'n! Cat/f as 
(—1)*n!/f /(n!/f). 

If we take the first of these to be 2, and the second 2, then 

> >> 
3, = and % = 
(a!) (n!/f) 


Then the square of the distance between C,! and Cyt); is (—1)"-*{n!}2/ (n!)? 
(n!/1)- 

Since this is independent of f, it follows that C,1 is the center of a circle pass- 
ing through the points Cain, Catz, » Calyn. It follows that C1234 passes through 
Cizs, Ciza, and is therefore the Morley circle of the 4-line; Ciss4s passes 
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through Cisse, Cizas, Cisas, and is therefore the Morley circle of the 
5-line; and so on indefinitely. 


We also note that the square of the radius of the circle C,1/z is 
{n!/f}? 
(n!/f)? 
3. The de Longchamps points.* Consider the circles Cis3 and Cig: 
(123)2 — z >. ky1*(23) — 2123 >> = 123 >> ke3(12) 
(124)23 — z >. k1*(24) — 2124 >> k,24(24) = 124 >> 
Eliminating 2, we have 
(123)z — 123 >> D> + 123 D> 
(124)z — 124 k,24(24) § + 124 ki keA(12) 
This reduces to 
[2(12) + 12{ — kel} >> — 1234 >> bi ke34(12)(34)] = 0. 


This means that these circles intersect at Ai: and also at a point which we 
shall call Kiss, which is given by 


1234 kiks34(12)(34)_ 


Evidently Cis, and C23, also pass through this point. This is the de Longchamps 
point of the 4-line. 


We now set up a chain of points K,1, defined as 
(—1)"n! >> kiken!/12(m!/12) — 229-6} 
> 


If we use for 2; and %, the values given in the last section for C,!, and for 
% the value given above for K,!, then 


(—1)* > kil (nt/1f)- 


> Ail (n!/1) 
Then 
21 
(n!/f) >> (n!/1) 


Hence the square of the distance between these two points is 
* de Longchamps, Joc. cit. Goormaghtigh calls these points the centric foci of the n-line. 


| 
| 
| 
| 
| 
| 
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(n!/f)? 
Comparing this with the result given at the end of section 2, we see that 
K, lies on the Morley circles +++, Thus Cress, 


C1348, Cosas, all pass through Kyi2345; and so on indefinitely. 
4. Clifford points and circles. | We shall not be concerned greatly with these. 
However, if we take the point Fizs«s to be 
—12345 k,k:345(12)(345) 
>> 119(2345) 
we may readily prove by the methods used above that Figs is the center of a 
circle passing through Kiess, Ki2«s, Kisas, This is the Clifford circle 
of the 5-line. The square of its radius is 
— {12345}2[ >> &1k:12(12)(345) }* 
[ |? 
It may readily be verified that this circle also passes through Kissss, which is 
Bath’s theorem. f 
Similarly it may be shown that Fizsise, which is given by 
— 123456 > kik2ks456(123) (456) 
>> 
lies on this circle and also on the Clifford circles Fiesss, Fissse, Fisase, 
This can be generalized to prove the existence of the Clifford circle of the 
2p+1-line and the Clifford point of the 2p-line. 
5. Miquel circles. We call the circle concentric with the Morley circle of the 


n-line and passing through its de Longchamps point, the Miquel circle§ of the 
n-line. Taking C,! to be 2: and K,! to be 2g, we find that 


21 — 22 

(n!) >> 


Hence the square of the radius of the Miquel circle is 
(a!)? 


t Clifford, Mess. Math., 1871. 

¢ Bath, Proc. Camb. Phil. Soc., 1939, p. 518. 

§ Miquel, Jour. de Math., 1838, p. 485, announced important properties of the 4- and 5-lines. 
He did not make explicit reference to these circles. 


| 
| 
| 
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Note that when  =4, the Miquel circle coincides with the Morley circle. 
6. Hermes points. Consider the circles Ci234 and Cis. 


(1234)2a—z k115(234)+21234 k:234(234) = — 1234 >> kik234(34) 
(123)22—2 k11(23) —2123 >> =123 
Eliminating 2, 
(1234)z + 1234 b:234(234) — 1234 hike34(34) {14 — 24} 
(123)z — 123 #:23(23) 12(23) + 123 {12 — 27} 
This reduces to 


[z >> k118(234), — 1234 >> ky k234(34)(12) 
[242(123) — 123 >, k:23(23) {12 + = 0. 


The Miquel circle of the 4-line thus intersects the Morley circle of the 
3-line obtained by omitting /, in the de Longchamps point of the 4-line and also 
in the point given by 

123 >| k:23(23){1? + 4°} 
4°(123) 


We shall call this one of the Hermes points of the 4-line and denote it by Hg; 12s. 
There are four Hermes points of the 4-line, H4:123, H3:124, H2:134, Ai: 234.* 
We shall now consider the point Hy: .1;; which we shall take to be 


(—1)"n!/f >> km!/1f(n!/1f) {12 + f?} 
f?(n!/f) 


We wish to show that this represents a Hermes point of the n-line, which we 
define as the intersection other than K,! of the Morley circle of the n—1-line 
€,1/7, and the Miquel circle of the -line. 

We take 2; to be the value just given. Then 


= > (n!/1f) {12+4?} /(n!/£). 
For zg: and 2, we take the values for C,./7 given in section 2. Then 
(—1)"{n!/f}? >> (n!/If) _ 
Then the square of the distance between Hy. »t;y and Cytyz is 
(—1)"{n!/t}? >> kil(m!/1f)- >> 
(n!/f)? 


* Hermes, Nouv. Ann., 1859, p. 359. 
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But this is the square of the radius of €,1!,;. Hence Hy, «1;; lies on this circle. 
Again, let 2, and % have the same values as above, but give 2: and % the 
values given in section 2 for C,1. Now 


21 — 22 + 1 — &s 


f?(n!) (n!) 
Then the square of the distance between Hy, »1;; and C,1 is 


(n!)? 


But this is the value obtained in section 5 for the square of the radius of 
the Miquel circle with center C,!. Hence Hy. atyy gives the » Hermes points of 
the n-line by taking f=1,2,---,m. 

7. The new chain of circles. We are now ready to set up the new chain of 
circles which it is the main purpose of this paper to announce. 

These circles have Hy: ,1;; for center and pass through Katy. 

We must first find the square of the radius of this circle. Here 2; and 2 have 
the values given in the last section for Hy, »1; and 


(—1)""'n!/f >> {12°-* — 229-8} 


22> 
> 
— kikel2(m!/12f) — 229-8} 
Ze = 
>> (n!/1f) 
Then 
and 
2. = 


(nt/f) kil 
Then the square of the radius is 
(—1)"{n!/f}* kil (m!/1f) {1? + £2} - { (1? + £7) } 
f°(n!/f)? 
We shall now show that this circle also passes through the points Hy, at/z,, 


where g is any value from 1 to m other than f. 
Here z; and 2, have the values used above, while 


(—1)""'n!/fg >> { 1? f?} 
f*(n!/fg) 


= 
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_ {12 + £2} 


(n!/fg) 
(—1)"{n!/fg}? >> kil (m!/1£) {12 + £2} 
f?(n!/f) 
g? { 12 +4 f?} 
(n!/f) 


And the square of the distance between Hj: »1;; and Hy. ,\yy, turns out to be 
the same as that given above. Hence we have proved 


THEOREM 1. Each of the n Hermes points of an n-line, such as Hy: n\/;, f being 
any integer from 1 to n (where the Miquel circle whose center is C,\ is intersected 
by the Morley circle Cy\);), ts the center of a circle which passes through the de Long- 
champs point K,\z and also through n—1 Hermes points Hy: Hy: 


’ 


Using these methods we may show that the equation of the circle with 
Ay: as center and passing through is 


(n!/f)f%3—2f? {1°+-£2} —2(—1)*n!/f >> km!/1f(n!/1f) {1°+£?} 
=(—1)"n!/f >> { 22-8} { 124-7} {22+ £2} 
In the same way, the circle with H,. ,1;, as center passing through K,1;, may 
be written down by replacing f by g in the above. 
Multiplying the first of these equations by g?(1f)(2f) - - - (fn) and the second 


by f?(1g)(2g) - - - (gn) and subtracting, we obtain the equation of the radical 
axis of these two circles to be 


of?g? >> — 2(—1)"{n!/fg}? >) 
= (—1)""'n!/fg >> kiken!/12fg(m!/12) — 22-8} {1424 — figs} 


Similarly, the radical axis of the circles of this system with centers Hy, a1/s 
and H;,; ,!;, can be written down by substituting h for g in the above expression. 
Eliminating 2 between these two equations, we find that 


_ kiken!/12(n!/12) {129-8 — 229-8} 
kil?*—"(n!/1) 
Since this is independent of f, g and h, it follows that 


THEOREM 2. The n circles of Theorem 1 form a system which has a common 
radical center. 


8. The 5-line. In the special case m =5, the radical center becomes 


312 } 
| 
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—12345 ki ks345(345) {12 — 2%} 
>> &:13(2345) 


But this is the value of Fiss4s given in section 4. 


THEOREM 3. The radical center of the five new circles of the 5-line is the center 
of the Clifford circle of the 5-line. 


Moreover it may be proved that 


THEOREM 4. The Clifford circle of the 5-line is orthogonal to each of the five new 
circles of the 5-line. 


It may also be proved laboriously that: 


THEOREM 5. The inverse points of the de Longchamps point of the 5-line with 
respect to each of its five new circles all lie on the Clifford circle of the 5-line. 


The equation of the straight line A12Cis, is easily shown to be 
+ 42} — kel {2* + 4°} ] 
(12) 


This may be proved to pass through 123. Similarly AisCisg and As3Cos4 pass 
through H4:123. This theorem of Hermes is the reason for giving the name of 
Hermes points to this set of points. 


4%; — 1°2% = 


COROLLARY. A12Ci23 passes through and Hy; 125. 
passes through H4:123 and Hg: 125. 
passes through Hs: 123 and Hy:124. 


There are three such lines passing through each of the ten vertices of the 
5-line. 
Again the equation of the straight line H4:123 Hs: 123 is 


1 k 2 42 12 52 
4°5%z + 12293%% — 23 + } 


similarly the equation of the straight line A; 12475: 124 is 
124 >> k:24(24) {12 + {12 + 57} 
(124) : 


These two lines intersect at —1234 ) hk, 234(234) {12452} /(1234), which is 
Hs, 1234- 


It may be shown that 134F2:134 and Hs: 234Fi:234 also pass through Hes, 134. 
This may be generalized as follows. 


+ 1°2°4°2 = 


THEOREM 6. Through Hy:syz pass four straight lines joining Hy. 517, and 


Hg: 5/40, where f has any of the values 1, 2, 3, 4 or 5, and g any of these values other 
than that taken for f. 


} 
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CoROLLARY. Hy: st)7g and H,; 51;7q lie on the line of centers of the new circles with 
centers Hy, and H,: 51/9. 
9. The 6-line. In the special case »=6, the radical center becomes 
123456 k1k:3456(3456) {14 — 24} 
>> k115(23456) 


It is easily proved by the methods used above that this point lies on each of 
the six new circles of the 6-line and also on the Morley circle of the 6-line. 


THEOREM 7. The six new circles of the 6-line intersect at a common point, which 
also lies on the Morley circle of the 6-line. 


Again, Ai; 1345 is 


345 k45(45) {12 + 32} {22 + 32 
1°22z + 374°5°g = v45(45) {1* + 3°} { +34} 


(345) 
and 1346 is 


(346) 


These intersect at 
—3456 >> k456(456) {12 + 32} {2? + 3°} 
1°22(3456) 


Evidently FN: 235¢L22: 1356 and Ai: 2as6ft2: 1486 also pass through this point. But this 
point may also be proved to lie on the new circles with Hii-23456 and He:13456 as 
centers. Hence 


THEOREM 8. In a 6-line the two new circles with centers Hj:23456 and Ho: 3456 
meet in a point other than that common to all six circles and through this point pass 
also the four lines Hy: 234sHs: 1345, 1346, Hi: 1386, Ai: 1486. 


This may also be generalized. 


It was proved in my previous article* that A12Kio4s, A1sKis4s, AosKosas inter- 
sect at a point which we shall call X123,4s, which is on the Clifford circle of the 
5-line. There are ten such points on this circle. 

Similarly, AisKisas, intersect at and Ai2Kisse, 
A13Kisse, intersect at X123,56. 

Now in the Clifford circle with center Fi2343 which passes through F, the 
Clifford point of the 6-line, 


Z X 123,48 = Z 123,45 = Z 


And in the Clifford circle with center Fiesae, 


* Loc. cit. 
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Z X 123,46 = Z 128,46 = Z 12. 


But Au, Ky234, Kia all lie on the circle Hence Z 
= Z Ki234FX 123,46. Therefore X123,45X 123,46 passes through F. Similarly X123,5s may 
be proved to lie on this line. We thus obtain 


THEOREM 9. The sixty points on the 6 Clifford circles of the 5-lines obtained by 
omitting in turn each of the sides of the 6-line, lie three by three on twenty straight 
lines all of which pass through the Clifford point of the 6-line. 


The analytic proof of this theorem is very laborious. Lachlan, in the preface 
to his Modern Pure Geometry* said: “Although Analysis may be more powerful 
as an instrument of research, .. . it might well be taken as an axiom, based 
upon experience, that every geometrical theorem admits of a simple and direct 
proof—by the principles of Pure Geometry.” If the writer had possessed enough 
acumen to develop his other theorems by such methods, a good deal of space 
might have been saved. 


* London, 1893. 


THE SOLUTION OF A SPECIAL HEAT AND DIFFUSION EQUATION 
DOUGLAS ALLAN, University of Toronto 


1. Introduction. The solution of the problem discussed here has been used 
to obtain results of some importance to geophysics [1-3]. Although the deriva- 
tion has little mathematical novelty, it has never been published, and may prove 
interesting, if only as a slightly out-of-the-ordinary example for differential 
equations courses. 

The equation to be considered arises in connection with the problem of 
finding temperatures and heat flows due to radioactivity within a sphere. It 
is shown in [2] that, assuming spherical symmetry, constancy of the thermal 
conductivity, etc., the required solution may be obtained by the superposition 
of a number of solutions of the same general type, namely, solutions of the 
equations 


with boundary conditions 
O(a, t) = 0, 
(2) 6(0, t) < @, 
6(r,0) = 0. 
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Call the temperatures in the regions 0SrSr, and Sra @, and @ respec- 
tively, and let 


(3) “= 
(4) = 70>. 
Write k for K/pc and J for Ho/pc. Then the problem is to solve the equations 
Ou 3? 
(5) 
ot or? 
(6) —=k + Jre™ 
ot or? 


subject to the conditions 
t) = 0, 0) 0, 


) = — 
—— 41\71, = — 
or 1\71 ar 2\"1 


u:(a, t) 0, u2(r, 0) = 0. 


The most obvious method of solution is to apply the Laplace Transform, 
and this will be given here. A second method, which gives the final solution in a 
slightly neater form, will be discussed a little later. 

2. Solution using the Laplace transform.* The subsidiary equations obtained 
by applying the Laplace transformation to (5), (6), and (7) are 


(8) 
dr? 
Jr 
(9) pie =k 
where 
#,(0)=0, =0, 
(10) fii(ri) = 


dr “dp 


The general solution of (8) is 


(11) = + 


* The notation used is the same as that in Carslaw and Jaeger: Conduction of Heat in Solids, 
Oxford, 1947. 


| 
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where g stands for ./p/k. The boundary condition @,(0) =0 will be satisfied if 
this is taken in the form 


(12) a@, = A sinh gr. 
Similarly, the general solution of (9) is 
rJ 
+) 


In order for the remaining boundary conditions (10) to be satisfied, the following 
equations must hold: 


(13) fl, = Bsinh gr + C cosh gr + 


aJ 
(14) B sinh ga + C cosh ga + ————- = 0 
+ A) 
(15) A sinh 
sinh gr cosh gr; + ————- = A sin : 
= 
(16) B h gr; + Cq sinh gr; + ~ A h 
cosh gr sinh gr, + —————- = Aq cos » 
q q +d) q 
Hence, 
J sinh _ 
om - — r, cosh g(a — — 
+ A) sinh ga q 
J cosh ga sinh gr; 
(17) B=- — cosh ga cosh gr; + ) 
p(p + A) sinh ga 
J sinh 
C= ( r, cosh an). 
b(p + A) q 


Thus the solutions of (8) and (9) satisfying (10) are 


J sinh sinh g(a — 
(18) a4, = + cosh g(a — + 
+ A) sinh ga q 
J sinh gr sinh gr; cosh ga 
i, = (« — r, cosh ga cosh gr; + ) 
P(p + d) sinh ga q 


(19) 


sinh gr: 
+ cosh gr — 4; cosh ar). 
q 


The solutions for u; and u, may now be obtained by inverting the above 
solutions, using the Fourier-Mellin theorem. The calculation of residues in- 
volved is tedious, but elementary, and need not be given here. The poles of the 
functions are at the points p=0, p= —A, and p= —m'x*k/a?, form=0,1,2,---. 
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Thus the solution obtained is in each case in the form of a single first term plus 
an infinite series, viz.: 


sin nr (1 — — cos n(@ — 71) — — sin (a — 11) 
(20) na a n@ 
(—1)™ — —cos 
2a? 2 + mar a a mr a 
(21) — =———— sin nr — —sinna+ (= cos — — sin nr; } sin (a — 1) 
Xsin na a a na 
(—1)" — — cos —— + — sin 
2a? a a mr a 
De sin ——. 
TR ad — 


where 7 stands for »//k. 
The solution for @ is therefore* 


sin — — cos — — —sinn(a — 
a na 
0 ae™ for0OSrsnj 
sin nr — — sin na + | — cos nr; — — sin nr; } sin n(a — 1) 
a a na 
(22) 
form, Srsa 
(—1)" — —cos + —sin 
2 a a mr a 
> (m' a") t. sin) ——. 
TRE a m(a?n? — 


3. Second method of solution. The second method of solution} is simply 
to try a solution of equations (5) and (6) of the form 


(23) = Ri(r)e™, 
(24) ui = R2(r)e. 


* This solution was first obtained by G. Comenetz [1]. That for the simpler problem with 
r,=0 has been worked out independently by H. C. Urey, and is given in his book, The Planets, 
p. 52. Both can be obtained quite easily from a very general solution by A. N. Lowan, Phys. Rev., 
vol. 44, 1933, p. 769. 

t This was kindly suggested to me by Professor G. F. Duff. 


; 

| 
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Then R; and R; satisfy 


(25) 
dr? : 
(26) + 1°R: + Jr = 0, 
dr* 
where 
R,(0) = 0, Ri(a) = 0, 
Ri(ri) = R2(n1), 
(27) 


dr 


These equations can be solved in much the same manner as equations (8) and 
(9), giving 


Jasin nr 
(28) Ri = — 1 — — cos n(a — r;) — — sin n(a — 14) 
A sin na a na 
r 
(29) R2 = — sin nr — —sin na + | — cos nr; — — sin nr; } sin n(a —7) }. 
A sin na a a na 


The corresponding solutions u{ and uj satisfy (5) and (6) and all of the boundary 
conditions (7) save the two 


u(r, 0) = 0, 0) = 0. 


However, this is easily corrected, for if 


(30) u=u'+ U, 

where 

(31) Osrsn 
uz nsrsa, 


then in order for u to satisfy (5), (6), and (7) completely, U must satisfy 


(32) —=k—», 
ot or? 
where 
U(a, t) = 0, U(0, t) = 0, 
(33) (a, (0, #) 


U(r, 0) = — u(r, 0). 


This of course is a standard heat conduction problem, and the solution is 


| = 
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(34) U= Ane sin =, 
m=1 a 
where 
2 mar 
(35) An = ~f U(r, 0) sin dr. 
avo 


The calculation of the A,, is fairly easy, and the final solution obtained is 
the same as (22), where now the summation term corresponds to U and the 
first term to u’. 

However, this method of solution points to a neater and better way of 
writing (22), for from (33) 


u'(r, 0) = — U(r, 0). 
But 
u'(r, t) = u’(r, O)e* 
= —U(r, O)e™ | 
(36) marr 1 | 
mar | (—1)" — —cos + —sin 
= > e*-sin . a a mr a 


m(a?n? 


Hence the final solution can be written 


> we marr 
= = — sin 
J ST ami a 
37 r 1 mrr 
(37) (—1)" — —cos 
. a a mr 
m(a*\ — m*x*k) J 


(remembering 
4. Remarks. The heat flow can at once be found from (37), viz. 


00 
F(r, 1) = —-K— 
or 
2a*Hok 1 
(38) = (e* —_ x" k/a*)?) . (-. sin — cos 
es r a a a 
1 mrr 
a mr a 


m(a*\ — m*x*k) 


The average temperature is also of interest. This is given by 


| 
| 
| 
| 
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3 
§= — 


4ra*® 0 
2 
m=1 
r mrr 1 
a mr a 


— m*x*k) 


The behaviour of these solutions is quite complicated. Some results for the 
special case of the Earth are given in the references; [3] contains a discussion of 
the reasons for using this solution, and the method of applying it to fairly com- 
plicated distributions of radioactivity. 

5. Application to the diffusion of argon. An interesting application of the 
above solution is to the problem of the diffusion of a gas being produced within 
a solid by radioactive decay.* The differential equation which applies is (using 
the notation for the case of the gas argon-40 being produced by the decay of 
potassium-40) 


where A* denotes the concentration of argon-40, Kq° the initial concentration 
of potassium-40, and \ and ), are constants. 

The solution required is that for the average argon content at time ¢ (denoted 
by A*°) of a spherical solid with the potassium distributed homogeneously 
throughout, and with zero argon content at t=0 and zero surface concentration 
of argon (i.e., similar boundary conditions to those of section 1). This solution 
is given by setting r,; =0 in (39), and substituting A“ for 8, D for k, and \Kq for 
J. Thus, 


A“ 6a? ert — 
m*(m*Dx*? — a*h) 


or, writing K* for c for and s for ax/X/D, 


(41) 


( ) c K* 
m? m? 1 
s? 


The diagram shows the behaviour of this function for various values of s. (The 


* This was pointed out to me by Dr. H. A. Shillibeer, who developed the results independently 
in connection with his researches on the potassium-argon method of dating minerals. The diagram 
and most of this section are taken from [4], with his kind permission. 


| 
| 
} 
i 
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X used is the decay constant for potassium-40, .55X10-* years!) Evidently 
there are two distinct types of curves. This can be seen from the formula (42), 
for 


t 1 ifs< 
43 lim { — )- m=1 mT $s 
( ) im (— K* s? 
+ 2 ifs > x. 


The limiting case s =7 is given by 


1 6X 
(44) — —=— 


dead 
c K* — 1) 


which, as the graph shows, is asymptotic from below to a straight line, viz. 
(45) 1 A® 
c K® 


It may be noted that if s were <7 in physical fact, it would be almost im- 
possible to use measurements of A‘*°/K*® to determine #, the age. However, 
experimental results obtained for argon seem to show s is large for most mineral 
samples, so that this age-determination method is practicable. 
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MATHEMATICAL NOTES 


EpiTep By F. A. FIcKen, University of Tennessee 


Because of the large number of papers on hand, consideration of new papers for this 
department has been temporarily suspended. 


EVALUATION OF PRODUCT SEQUENCES BY MATRIX METHODS* 
GEORGE BRAUER, University of Minnesota 


1. Introduction. Let 


(1) on = 
k=O 
be a regular transformation by means of which a sequence So, si, - - - , is evalu- 


able to a if the series in (1) converges for each m and ¢,—0 as n— ©. The hypoth- 
esis that (1) is regular means that ¢,—s whenever s,—s. We shall use the well 
known conditions 


(2) | SM n= 0,1, 
k=O 
(3) lim = 0 k= 0,1,2,--> 
(4) lim > ane = 1, 
k= 


which are necessary and sufficient for regularity of (1). An exposition of these 
matters is given, for example, by Hardy [1]. 

We shall say that the transformation (1) and the matrix A =(a,,) have the 
product property P if A is real, regular, and such that if s, and s, are real 
bounded sequences evaluable to o and o’ respectively, then the product sequence 
Sn$x is evaluable to the product ao’. 


* A similar problem is considered by Mazur and Orlicz: On linear methods of summation, 
Studia Math., vol. 14, 1954, pp. 129-160. 


| 
| 

| 
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2. Sequences of zeros and ones. The following theorem shows that the 
methods of Cesaro, Abel, and Borel and many other standard methods for 
evaluation of sequences fail to have property P. Moreover it shows that if A has 
property P, then relatively few sequences of zeros and ones are evaluable A. 


THEOREM 1. If A has property P and if so, S51, > + + , is @ sequence of zeros and 
ones evaluable to o, then o is 0 or 1. 


To prove this theorem, we note that if A has property P and if s, is a 
sequence of zeros and ones evaluable A to a, then the sequence s,s, (which is 
the same as the sequence s,) is evaluable A to o?, so that o=o? and a is 0 or 1. 

We now employ some results of which an exposition has been given by Hill 
[2]. To each sequence So, 1, S2, - - + of zeros and ones there corresponds a point 
t of the interval 0 S$¢<1 having the binary expansion sos:52 - - - . If A is regular, 
then the set of sequences of zeros and ones evaluable A corresponds to a set of 
values of t having Lebesgue measure 0 or 1, and, moreover, if the measure is 1 
then the set of zeros and ones evaluable A to 1/2 corresponds to a set of values 
of t having measure 1. It follows from this and Theorem 1 that if A has property 
P, then the set of sequences of zeros and ones evaluable A corresponds to a set of 
values of ¢ having measure 0. 

3. A criterion for property P. 


THEOREM 2. A real regular matrix A has property P if and only if s2 is evalu- 
able A to o* whenever s, is a real bounded sequence evaluable A toc. 


Let a real regular matrix have property Q if s? is evaluable A to o? whenever 
Sn is a real bounded sequence evaluable A to o. Since we may take s,/ =Sa, it is 
obvious that if A has property P then it also has property Q. To prove the con- 
verse, let A have property Q and let s, and s, be real bounded sequences 
evaluable A to o and o’ respectively. The real bounded sequences {s,+s, } 
and {s,—s, } are then evaluable A to ¢+o’ and o—o’ and hence the sequences 
{(sa+se)*} and {(sa—sy)?} are evaluable A to (¢+0’)? and (¢—o’)?. Since 
AsaSn )? and 400’ it follows that 
is evaluable A to go’ and Theorem 2 is proved. A variant of Theorem 2 which 
may be useful is the following. 


THEOREM 3. A real regular matrix A has property P if and only if s2 is evalu- 
able A to 0 whenever s, is evaluable A to 0. 


This is a consequence of the fact that if s? is evaluable A to 0 whenever s, 
is evaluable A to 0 and if s,’ is evaluable A to o’ then {s,’ —c} is evaluable A 
to 0, so {(s,, —c)*} is evaluable A to 0, and linearity of A implies that (s,/)? is 
evaluable A to (¢’)?*. 

4. Consequences of property P. 


THEOREM 4. If a real bounded sequence s, is evaluable to o by a matrix having 
property P, then 


| 
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(5) lim inf | |o| limsup | 

To prove Theorem 4, it will be sufficient to prove that 
(6) lim sup =| 2 |e — x2 0, 


because putting x=0 in (6) gives the last inequality in (5), and choosing x so 
great that x >o? and x>s? for each n gives the first inequality in (5). To prove 
(6), suppose that there is some real x for which it is false. Then xo? and there 
is a number @ for which 


— ~ 
(7) lim sup <0<1. 


With M satisfying (2), choose a positive integer m such that 6"M <1. Let 


2 
Sa & 


(8) t, = 


o— 


The hypothesis of Theorem 4 implies that s2 is evaluable A to ¢? and hence that 
s?—x is evaluable A to ¢?—x so that ¢, is evaluable A to 1 and & is evaluable A 
to 1. Thus 


(9) lim = 1. 


Because of (7) and (8), we can choose ky such that | ts | <6 when k2kp. Use of 
(9) and (3) then gives 


(10) lim = 1. 
keke 
But use of (2) gives, for each m=0,1,2,+--, 


= m 
> Onkt 


(11) 


1 
S06 D | ans | 


The contradiction in (10) and (11) establishes (6) and completes the proof of 
Theorem 4. 


THEOREM 5. If a real bounded sequence s, is evaluable to o by a matrix A having 
property P, then o is a cluster value (or limit point) of the sequence S,. 


The hypothesis implies that the sequence so—o, s1—¢, - + + is evaluable A 
to 0 and hence Theorem 4 implies that 


lim inf |s, —¢| SOS limsup |s, —o|. 


ae 
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The result follows. 
5. Conclusion. For the following remarks, the author is indebted to R. P. 
Agnew. The transformation defined by o9=50, and 


{ Tan = — + n> 0, 


= — 45en—1 + n> 0, 


(12) 


has property P. This is an easy consequence of the fact that ¢,—>L if and only if 
San—L and S3,-1—L, the values of the numbers $3,-2 being completely irrelevant. 
It may be possible to show that each transformation having property P must 
have a form resembling (12) in essential respects. 

The following seems to be an entirely reasonable conjecture. A given matrix 
and transformation A have property P if and only if there is an increasing se- 
quence ™, m2, m3,°-- of non-negative integers which depends only upon A 
and which is such that an arbitrary sequence So, si, S2, -- + is evaluable to L if 
and only if the special subsequence Sa2, Sas, * Converges to L. 
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RATIONAL UNIVALENT FUNCTIONS 
W. C. Royster, Alabama Polytechnic Institute 


Let be the class of functions f(z) =s+ 9-2 which are regular and 
univalent in the unit circle E(|z| <1). The subclasses of ¥ whose coefficients a, 
belong to a quadratic field have been studied by Friedman [1] and Bernardi [2]. 
Recently Linis [3] gave a short proof of a theorem of Friedman which states 
that if the a, are rational integers then f(z) is a rational function and is one of the 
nine functions: 


z,2/(1 +2), 2/(1+ 2%), 2/1 +2)*%, 2/1 +2+2°). 


He also gave some rational functions whose expansion had Gaussian integers as 
coefficients. 

The purpose of this note is to extend the method of proof given by Linis to 
quadratic fields with negative discriminant and to give criteria for selecting 
which combinations of possible coefficients yield rational functions which are 
regular and univalent in E. 

For completeness we will sketch Linis’ proof. Let a, be a rational integer. Let 
2/f(z) = b,2". By Prawitz’s theorem [4] (n—1) | Clearly 
bo =1, = and since f(z) is univalent in E, | so that bi | <2. From 
Prawitz’s theorem it follows that |b,| $1, 22. The coefficients }, can be com- 
puted by the determinant (—1)"b,=A(C,;) where Cy=0, j7>i+1; Cij=1, 
j=tt1; 4, 7=1,---, Hence the are rational integers. The 


= 
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possible values of b, are: bb =1; b:=0, +1, +2; b2=0, +1; 5,=0, n23. Certain 
combinations of these values yield the nine above-mentioned functions. 

Let R(./m) be a quadratic field where m is a square-free rational integer. 
If m=2 or 3 (mod 4) the integers of R(./m) are given by a+b\/m, where a 
and b range over all rational integers. If m=1 (mod 4) the integers of R(./m) 
are given by }(a+b./m) where a and b are both even or both odd rational 
integers. We will be interested only in the cases in which m is one and in which 
m is negative. The units of R(./m) are those numbers whose norm is +1. 
The only units of R(./m) where m is negative are +1 unless m=—1, —3, 
in which case the additional units are +i and +w, +a, respectively, where 
w=4(—1++/31) and @ is the complex conjugate of w. (For the arithmetic of 
algebraic number fields see [5], Chap. VII). 

Let a,€ R(/m), for a given m, the b, are given by the previously mentioned 
determinant, so that },€(Rv/m). Since R(./m) has no elements whose modulus 
is a number between zero and one, it follows from Prawitz’s theorem that d,, 
n= 3, in the expansion of z/f(z) are all zero. 

The possible values of 5, are as follows: 


m=—1; bo=1; 6 =0, +1, +2, +i, +24, +(1 9); 
= 0, +1, +3. 

m=—2; bo=1; 6 =0,+1, +2, +72i,1 + 24, 
bs = 0, +1. 


m=—3; b=1; 0b =0,+1, +2, to, +0, +20, +20, +(1 
+(1 — &), — &);b, = 0, +1, tw, ta. 


m< —3; bb =1, =0, +1, +2; b:=0, +1, the same as when m=1. 
Before enumerating the }, which yield regular univalent rational functions, 
let us state two lemmas. 


Lemma 1. The function f(z) =2/(1+b:2) ts regular and univalent in E if and 
only if | bs st. 


LemMA 2. The function f(z) =2/(1+bie+bse*) with =1 is regular and 
univalent in E if and only if b,/(b2)'!? is real and in modulus less than or equal 
to 2. 

The proof of Lemma 1 is clear. The univalency in Lemma 2 follows immedi- 
ately. Since |6:| =1, f(z) will be regular in E if and only if the zeros of the 
denominator lie on the periphery of E. Let z=s/(b:)'/?, then the denominator 
is s*+ [bi /(b2)"/?]s+1. Its zeros are on | s| =1 and hence are complex conjugates. 
Using these properties and the quadratic formula we see that the zeros will be of 
unit modulus if and only if },?/b: is positive and less than or equal to 4 which 
implies that b,/(b:)'/? is real and in modulus less than or equal to 2. 

We note that every field R(./m) includes the field R(./1), hence we will 
enumerate those 5, not concerned with R(+/1). 


| 
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Making use of the following facts: (w)?=—a, (—w)?=ia, (a)/?=w, 
(—@)/? = —iw, and the properties of the cube roots of unity, we can easily con- 
struct a table of values of the 5, which satisfy the two lemmas. 


bs 0 1 -1 -i 
bi, m= —1 +i none | +4, | 0, | 0, (1-1) 

bi, m= —2 none none + V2i none none 

b,m=—3 | +o,+@ | none| +(w—a) none none ta, +20) 0, t(1—w) |0, tw, +2] 0, + (1—a) 
bh, m<—3 Same as when m=1. 


This proves the following theorem first given by Bernardi [2] by other 
methods. 


THEOREM. The subclass of § having algebraic integral coefficients in R(./m), 
m=1 or a square-free rational negative integer, consists of forty five rational func- 
tions of the form f(z) =2/(1+b:2+-b.2"), where b, and bz are given by the above table. 
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A CONJECTURE CONCERNING NEIGHBORING TETRAHEDRA 
F. BAGEMIBL, The Institute for Advanced Study 


Two bodies in space are said to be neighboring if they have no interior 
points in common and there is a surface belonging to both their frontiers. A 
(finite or infinite) number of bodies are said to be neighboring if every pair of 
them are neighboring. 

It is known that it is possible to have infinitely many neighboring bodies— 
even infinitely many neighboring convex bodies.* This suggests the problem of 
defining, by means of some property, a class of bodies such that the maximum 
possible number of neighboring bodies of this class is finite, and of determining 
this maximum number. 

Consider the class of tetrahedra. It is not difficult to verify that at most 4 
neighboring tetrahedra can stand with their bases on one and the same plane 
(and lie, with the exception of these bases, on one and the same side of this 
plane). Consequently, the maximum possible number of neighboring tetrahedra 


* For the literature on this result, as well as an interesting discussion of some related prob- 
lems, see Heinrich Tietze, Geliste und ungeliste mathematische Probleme aus alter und neuer 
Zeit, Munich, Biederstein, 1949. 
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is certainly not greater than 17, and we conjecture that it is 8. That it is at least 
8 is shown by the following example. 


Cc 


A a B 


In the accompanying figure, all the lines are supposed to be coplanar. Let 
Q and g be points on opposite sides of the plane ABC. Then the following 8 
tetrahedra are obviously neighboring: 


QABE, QBCF, QCAD, QDEF, 
qabd, qbce, qcaf, qdef. 
ON ANALYTIC FUNCTIONS HAVING A POSITIVE REAL PART 
IN THE UNIT CIRCLE 
V. F. Cowtinc, The University of Kentucky 


A theorem of basic importance in the theory of univalent functions is the 
following: 


THEOREM. Let f(s) =1+ Doe, a,2" be regular and have Re(f(z))*20 in 
|z| <1. Then |a,| <2, n=1,2,---. 


* By Re(f{(z)) we understand the real part of f(s). 


JAIN 
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The function f(z) = (1+2)(1—z)—' shows that this bound on the coefficients is 
sharp. 

Proofs of the above theorem have been given by Carathéodory [1], Toeplitz 
[4], Nehari [3] and others. It is the object of this note to give another proof 
which is elementary in character and would seem to suggest from its structure 
the theorem to be proven. For this purpose we state the following lemma whose 
proof depends only on simple manipulations of the Cauchy integral formula. 
Such a proof may be found in Goursat [2]. 


Lemma. Let f(z) = U(x, y)+éV(x, y) be regular in | SR. Let z=re® and 
¢=Re*. Then for |z| <R we may write 


(1) fle) = iv(0,0) +~ 
0 


The proof of the above theorem now proceeds as follows. We have 


dy. 


fe) =1+ Saws, 


where a,=f™(0)/n!. But by (1) above 


f(m)(0) 
=— | 


n! 


where we take R<1. Also since U(0, 0) =1 the mean value theorem for harmonic 
functions yields 


1 
U(Re*¥) dy = 1. 
2rJ 


But since U(x, y) 20 it follows that 
1 1 Qn 
R*\a,| s f | U(Re*¥) | an = =2-1=2 
2rJo 2rJo 


for all R<1 and therefore S2,#2=1,2,---. 
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ON THE nTH DERIVATIVE OF COMPOSITE FUNCTIONS 
MicHEL McKIERNAN, Illinois Institute of Technology 
Let f be the function whose value at x is f(x). Let D*f denote the mth deriva- 


tive of f, that is, the function assuming for any x the value D*f(x). If f and g are 
two functions, then 


(1,) D(f +g) = D¥+D*g 

and according to Leibniz 
n| 

ri(n = r)! 


Let f and g be functions. Throughout this note the juxtaposition of two 
functions will denote the substitution of the second into the first (hence, the dot 
will not be omitted in products of functions). For example, (fg) (x) = (fg(x)) and 
(Df)g is the result of substituting g into Df. The derivative of the function 
fg is then given by: 


(31) D(fg) = Df)g-Dg. 


Formulae for D*(fg) have been found by Faa di Bruno and Schlémilch.t The 
purpose of this note is the establishment of two other formulae for D*(fg). The 
first is 


Formula (3,) gives an explicit form for the operators Y,,, in the relation 
(4n) D(fg) = Ya.r(g) 
r=1 


which was used by Schlémilch. It is known that the Y,,, satisfy the recursion 
formula:§ 


(S) Yotie(g) = Dg- + 
where 
(6) Yn,o(g) = O = Yuings(g). 


Let f be the function assuming at x the value f(x) =x*/k!, where OSkSn. 
Substituting this f in (4,) yields a second recursion formula 


t See, for example, Charles Jordan, Calculus of Finite Differences, p. 33, Sopron, Hungary, 
1939. 

§ Riordan, Derivatives of composite functions, Bull. Am. Math. Soc., vol. 52, pp. 664-667, 
1946. Also C. Jordan. 
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(7a. = — DY) — 
n,k n, g (k—»)! n,r\Q)- 
To prove (3,), one must show that 
r —1 re 
(8n,2) = >> g”*-D*(g*). 


=o SI(r—s)! 


The proof is by induction. By (3), Yia(g)=Dg. By induction, in view of (6) 
and (5), =D"g, which agrees with (8,1). Assume (8,,) for k=0,1,---, 
r—1. To prove (8,,,), in view of (7,,-), one has to show 


r —{)r- 1 rl r—k k —1)*- 
= — De) ae 


=o Si(r—s)! 
The term for s=r on the left equals D*(g’)/r! on the right. For a fixed s, 
0ssSr-—1, the coefficient of g*-*- D*(g*) on the right is 


r—1 (—1)*- 


bee (7 — — 5)! 


Let k=s+, then this becomes 
4 1)? 
pmo (7 —s— p)ipl 


which equals 


si\(r — s)! 


— — [sr — = 


But this is also the coefficient of g*-*- D*(g*) on the left. This proves (8,,,-). 

The second formula for D*(fg) is 

1 n! 

D(fg) = 2: ———— -- -(Dg) 
r=1 r| pi! pr! 

where the second summation is taken over the integral solutions of #:+ - - - 

+p,=n and p;>0, fori=1, - - -,r. The proof of (9,) is based on the following. 


LemMA: Y,,,-(g) consists of those terms in the expansion of D*(g*)/r! which do 
not contain g as factor. 


By (3), Yi(g) =Dg. Assume that Y,,-(g), r=0, 1, - - - , # contains no terms 
with g as factor. It follows by (5) that Y.4:,-(g) contains no terms with g as 
factor. But by 
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Yn,r(g) = g- Dg) + page), 
ao Si(r — s)! rl 


Hence the first term on the right, containing g as factor, must cancel with the 
corresponding terms in the expansion of D*(g’)/r!. 

The proof of (9,) is completed by applying the multinomial analogue of 
(2,) to D*(g’). 


1 1 n| 
r| pil: 


pe! 
where the summation extends over the non-negative integral solutions of 
fit +--+ +p,=n. A factor g appears whenever one of the derivatives is of order 


0. The stipulation p;>0 eliminates these terms, yielding, along with the lemma, 
formula (9,). 
In conclusion, we note the following consequences of the above formulae. 
In view of (5) it follows by induction that Y,,.(g) = [Dg]* for any n 21. Hence, 
(8n,n) yields the following identity 
n —1 
(10) 


g”*-D*(g*) for any g. 
=o —s)! 


It may also be shown that 


(n — 1) 


= — Dg: [Dg]. 


Hence: 
2(n—1)!2, 


(11) D*- = @ +i)! 


g”-*-D**1(g*) for any g. 


Also, comparing (9,) with the corresponding formula of Faa di Bruno, we obtain 


the identity 
i (D71g)- -- + -(D™g) 


(12) =—) 


where the summation is over where the summation is over 
and and 
Qi + 2gn + gn = pi > 0. 


Many cther identities may be found by particular choices of g in (10), (11), and 
(12). 


= 
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EpiTep sy G. B. THomas, Massachusetts Institute of Technology 


Because of the large number of papers on hand, consideration of new papers for this 
department has been temporarily suspended. 


A FUNDAMENTAL THEOREM OF CALCULUS 
M. K. Fort, Jr., The University of Georgia 


The following theorem is frequently called the fundamental theorem of the 
differential calculus. 


THEOREM A. If F, and F, are functions which have the same derivative at each 
point of an interval I, then there exists a constant c such that F,(x) — F2(x) =c for 
each x in I. 


The proof of Theorem A follows easily from: 


THEOREM B. If F is a function whose derivative is 0 at each point of an interval 
I, then there exists a constant c such that F(x) =c for each x in I. 


The usual proof of Theorem B is based on the Mean Value Theorem, the 
proof of the Mean Value Theorem is based on Rolle’s Theorem, and proof of 
Rolle’s Theorem is based on the Maximum Value Theorem for continuous func- 
tions. Since a proof of this last theorem is often omitted from the first course in 
calculus, a complete proof of Theorem B is usually not given. Recently, Pro- 
fessor Harry D. Ruderman pointed out to the author that it would be desirable 
to have a direct proof of Theorem B which would avoid the use of the Mean 
Value Theorem. In this note we give such a proof. 

We start with the fact that every bounded, monotone sequence of real num- 
bers converges. This fact is stated as an axiom in many calculus texts. It is 
then easy to prove that if JiDJ2>DJ;) --: is a nested sequence of closed 
intervals which is obtained by the “bisection process”, then there exists a unique 
number p which is a member of each J,. Moreover, if x, is a member of J, for 
each n, then the sequence %, x2, x3, : - - converges to p. We now give a proof of 
Theorem B which is based on these facts about nested sequences of closed 
intervals. 

It is sufficient to prove that if a and 6 are points of J then F(a) = F(b). We 
may assume that a<b. 

If J=[u, v] is a closed interval which is contained in J, we define S(J) 
=| (F(u) — F(v))/(u—»)|. The following lemma is obvious and its proof is 
omitted. 


Lemma. If u<t<v, J’=[u, ¢], J’ =[t, v] and J=[u, v], then either S(J’) 
= S(J) or S(J”) 2 S(J). 


We make use of the Lemma and the “bisection process” to obtain a sequence 
Ji +--+ of closed intervals such that J:=[a, 6] and S(Jn41) 2 S(Ja) 
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for each n. There exists a unique point p which is a member of J, for each n. 
We now use the Lemma to obtain for each m a non-degenerate subinterval K, 
of J, such that one end point of K, is p, the other end point x, of K, is also 
an end point of J,, and S(K,)2S(J,). The sequence xi, x2, - , converges 
to p, and since F’(p)=0 we obtain lim,.. (F(x,) — F(p))/(x,a—p) =0. Thus 
lim,.«. S(K,) =0. Since S(K,) 2S(J,) 2 S(Ji1) 20 for each n, we see that S(J,) =0. 
This implies that F(a) = F(b), which concludes the proof of Theorem B. 


SYSTEMS OF LINEAR DIFFERENTIAL EQUATIONS* 
W. G. Leavitt, University of Nebraska 


Let {y;} (¢=1, - - +, m) be a set of independent variables, and suppose that 
a system of m linear homogeneous differential equations with constant coeffi- 
cients is to be solved. The usual approach to this problem is to introduce de- 
rivatives of the {y,} as auxiliary variables, reducing the system to a (larger) 
set of equations all of the first order. In the following, an alternative approach 
is given, which makes use of the “elementary divisor theory” for matrices of 
polynomials. 

Using matrix notation, the system of m equations may be written 


(1) T(D)Y = 0, 


where T(D) is an m by nm matrix whose elements are polynomials in the differen- 
tial operator D, and Y is a column vector whose components are the 4;. 

A natural approach to the solution of (1), analogous to that of the theory 
of a single linear equation, is to assume a solution of form 


(2) = (Ko+ Kit +--+ + 
where the K, are constant vectors and a is a scalar constant. Substituting in 
(1) and equating coefficients of x to zero, we obtain a series of equations 
T(a)K, = 0 
(3) T(a)K.-1 + sT’(a)K, = 0 
s(s — 1) 
2! 


T(a)Ko + T’(a)Ki + T"(a)K2 +--+ + T(a)K, = 0. 


It is clear from the first equation that a necessary condition for the existence 
of a solution is that T(a) must be singular; that is, the equation T(a)X =0 must 
be solvable. (If T(a) is square, for example, this means that a must be a root 
of the characteristic equation | T(x) | =0.) 


T(a)K,-2 + (s — + 


T’(a)K, = 0 


* Presented (in part) to the Mathematical Association of America, Nebraska Section, May 2, 
1953. 
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In some casesf it is convenient to solve these equations iteratively for 
K,, K,-1, - - +. However this procedure leads in general to certain difficulties 
such as arise, for example, when the first equation of (3) has several solutions 
not all of which make possible a solution of the second equation, and so on. Fur- 
thermore, it is not clear from this approach just what value of s should be used, 
nor is it evident that all solutions must necessarily be of form (2). 

To analyze this problem more carefully we take advantage of the isomor- 
phism* between the polynomials in D and the set of all polynomials in an in- 
determinate \. This makes possible the use of a well-known theorem which states 
that there exist square matrices P(D) and Q(D) each having unit determinant 
such that U=PTQ has form 


(4) U(D) = uy 


0) 
in which the rank h is the same as that of T(D) and the { us} (the “invariant 
factors”) are so arranged that mu, divides uz, u_ divides u3, and so on.§ The 
matrices P(D) and Q(D) may, in fact, be obtained as products of elementary 
matrices. Thus substitution in (1) of Y=QX and multiplication by P leads to 
the equation 

(5) UX =0. 


The first result is that if h<m there will be »—h solutions of (5) involving arbi- 
trary functions. These will be of form 
(6) X; = (i= h+i,--- , n) 


with y, arbitrary and 6; the vector having 1 in the ith position, zero elsewhere. 
Using the usual properties of the differential operator, it is easily verified that 
(1) has corresponding solutions 


Vi Qio 


(7) Yi = + + 


(G=k+1,---,n), 


t See Frazer, Duncan and Collar, Elementary Matrices, Cambridge University Press, 1952, 
p. 172. 

* For simplicity, we will henceforth use this isomorphism without explicit reference. For exam- 
ple, when we say that g(D) is a factor of f(D) we will mean that g(a) divides f(A). 

§ See, for example, C. C. MacDuffee, Vectors and Matrices, Carus Mathematical Monograph 
No. 7, Math. Assoc. of Am., 1943, p. 143. (The generalization to non-square matrices is obvious.) 
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where Q% indicates the ith column of the kth derivative of Q(A) evaluated at 

Now suppose that one of the invariant factors, say u;, has a factor (D—a)* 
(an “elementary divisor”). Then clearly (5) has solutions 


(8) X; = te, (j =0,1,-++,¢—1). 


Again using well-known properties of the differential, the corresponding solu- 
tions of (1) are 


(9) = + jtQia + de 


where QY is the ith column of the jth derivative of Q(A) evaluated at \=a. 
Note that this is a solution of type (2) with s=/j. 

Now it is clear from an examination of the form of U(D), (4), that the only 
solutions of (5) (other than those given in (6)) are of this type. Then since every 
solution of (1) corresponds to a solution of (3), it follows that the set of all solu- 
tions of type (9), together with (7), give all solutions of (1). Hence a knowledge 
of the elementary divisors of T(D), and of the transformation matrix Q(D) 
permits us to write down immediately the general solution of the system (1). 


ON FINDING DIAMETERS OF CONIC SECTIONS 
E. A. HEDBERG, University of South Carolina 


Nearly all current textbooks in analytic geometry give special formulas for 
diameters of conic sections in which the chords have a given slope m, listing 
each type of conic separately and giving only the standard forms, where the 
center or vertex is at the origin. That all these formulas are related seems, from 
this writer’s experience, not universally known. I state the general case as a 
theorem. 


THEOREM. Let f(x, y)=Ax*+Bxy+Cy?+Dx+Ey+F=0 be the general 
equation of a conic section. The equation of the diameter which bisects the chords of 
slope m may be found by differentiating f(x, y) =0 with respect to x and replacing 
dy/dx by m. 


The statement of the theorem implies the existence of chords with slope m. 
It will be pointed out in the proof of the theorem that the condition that 
chords exist is that A+Bm-+Cm’#0. Attention is also called to the fact that 
the theorem applies to the degenerate conics. 

One may establish the theorem by direct substitution of y=mx+k in 
f(x, y) =0, setting up the formulas for midpoints of the chords, eliminating the 
parameter & from these expressions, and comparing the result with the theorem 
stated above. However, a proof which avoids such manipulations would seem 
desirable. The following proof involves a minimum of actual computation. 
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The fact that the locus of the midpoints of the chords with a given slope m 
is a straight line may be verified without determining the actual locus. For if 
one substitutes mx +k for y in the expanded form of f(x, y), in order to determine 
the abscissas of the points of intersection of the line y=mx+k with the conic 
f(x, y) =0, it is at once apparent that k does not enter A+Bm-+(Cm’, the coeffi- 
cient of x, and that it enters the coefficient of x to the first degree only. If 
A+Bm-+Cm!*<0 the quadratic equation has two roots, x; and x2, and the ab- 
scissa, X =4(x1+4%2), of the midpoint of the chord involves k to the first degree 
only. The ordinate, Y, of the midpoint, found from Y=mX-+-8, is linear in k, 
hence it is plain that the equation obtained by eliminating k from the expres- 
sions for X and Y is linear in X and Y. This straight line will henceforth be re- 
ferred to as the diameter. 

The coefficients which appear in the equation of the diameter are rational 
functions of m and the coefficients of f(x, y). The only denominators which ap- 
pear are constant multiples of A+Bm-+Cm?, hence we may express the equa- 
tion of the diameter using polynomial coefficients by clearing of fractions. The 
method used in obtaining the equation is obviously independent of the type of 
conic under consideration, the only limitation being that A+Bm-+Cm? is not 
zero. 

If A+Bm-+Cm?=0 then there is at most one finite point of intersection of 
the line y=mx-+k with f(x, y) =0, whatever value of k is used. That is, in this 
case, there are no chords, in the ordinary sense, for the corresponding value of 
m. The parabola, for instance, has no chords parallel to its axis, the hyperbola 
none parallel to an asymptote. In the following, therefore, it will be assumed that 
A+Bm+Cm’? <0. 

We now find the equation of the diameter, using the expanded form of f(x, y) 
without making any explicit substitutions. If mx+ were substituted for y in 
f(x, y) =0, and the condition for equal roots applied, which is also the condition 
that y=mx-+k be tangent to f(x, y) =0, the result would be (d/dx) f(x, mx +k) =0. 
However, 


d 


yum z+k 


That is, the points (x, mx+k) at which the lines of the system y=mx+k are 
tangent to the conic f(x, y) =0 lie on the straight line 


a a 
+ y) = 0. 


Now for an ellipse there are always two real points of tangency. The diameter 
plainly passes through these points, hence the above process leads to the equa- 
tion of the diameter. But we have seen that the coefficients in the equation of 
the diameter are independent of the type of conic involved. Thus, in every case, 
the differentiation process leads to the equation of the diameter. 


| 
| 
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An application of the above process to finding the equations of the medians 
of a triangle may be of interest. The following represents the solution in the 
general case. Let Lh: =Aix+By+C,=0, Ls=Axt+By 
+C;=0 be the equations of the sides of a triangle. Let us assume that B; +0. 
To find the equation of the median upon the side Z;=0 write the equation of 
the diameter of the degenerate conic L:L:=0 which bisects the chords of slope 
—A;/B;. Applying the general rule of the theorem above, one arrives at the 


— + L2(AiBs — AzB,) = 0. 


or 


In determinant form this equation becomes: 


0 
A; A: A; = 0. 
B, Ba Bs 


It may be readily verified that the last equation represents the median upon the 
side L;=0 even when B,=0. 


A COMMENT ON SNYDER’S CORRECTION FORMULA* 


EVELYN FRANK, University of Illinois 


The correction formula (2) of [5] goes back to Lagrange [4]. It is also found, 
for example, in [6], pp. 222-223, and has recently been repeated a number of 
times, for example, in [1] and [2]. In [3] it is extended to m equations in n 
unknowns. 
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ON THE APPLICATIONS OF THE FUNDAMENTAL THEOREM OF INTEGRAL 
CALCULUS 


Paut ScHILLO, University of Buffalo 


In the usual elementary calculus course, it is an exceptional student who can 
avoid the misconception that the definite integral is the sum of infinitely many 
infinitesimal differentials. This is because, in most texts, none of the problems 
requiring the application of the Fundamental Theorem of the Integral Calculus 
discourage this notion. In fact, the integrals, which yield the correct solutions 
to these problems, are most easily set up by using this misconception as a device. 
For example, the area of the region bounded by the line y=x and by the parab- 
ola y=x? can be thought of as the sum of the areas of all infinitesimal x —<x? 
by dx rectangles from x =0 to x = 1 and, consequently, can be written as: 


fe — x*)dx = }. 


Some texts attempt to remedy this situation by including among their 
exercises questions which ask for the values of limits like: 


1 n 
(1) Lim — >> and 
1" 
kr 
(2) Lim — ¥* sin — - 
no n 


Of all those who have, at some time, been introduced to the integral calculus, 
few would feel competent to evaluate these limits. Some might be successful 
after much effort in showing the validity of the following not very general for- 
mulas: 


n 2 1)? n k 
k=1 4 k=l n 2n 


Then, applying certain theorems on limits, they could show that 


12 2 
Lim — >> = rin 


1” 4n* 4 


and 


2n 
Lim sin = 2] Lim cos | Lim ———— | = 2. 


no n neo n 


However, the student who has not been weaned from the correct definition 


| 
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of the definite integral as the limit of a sequence of sums should be able to recog- 
nize (1) and (2) as the respective integrals: 


1 1 
f dx = — and f sin xdx = 2, 
0 4 0 


Thus, these limit problems help the student to learn just what the definite 
integral is. 

Now, according to a widely held belief, the methods used in solving the so- 
called “practical” problems (e.g., finding the area between curves efc.) are re- 
tained by students much longer than those used in solving straightforward 
drill problems of the same complexity. If this is true, much of the value of the 
limit problems is likely to be lost unless they are clothed in “practicality.” 
For this reason, an elementary calculus course might advantageously include 
problems of the “ultimate-average” type. 

Let S be a structure or situation, which can involve any arbitrary number 
of members of a certain kind; let A, be some arithmetic average associated with 
S when it involves exactly » members; then, if Lim,... A, exists, this limit is 
called an ultimate average in S. For example, the ultimate average distance 
between pairs of vertices in a regular polygon circumscribed by a unit circle is 
4/mr units, because 


and 


2 n ei. 4 
Lim A, = —| Lim Lim = 5) sin =]= —. 
avo n— 1 n 
Also, the ultimate average distance between pairs of points in a one unit column 
of equi-spaced points is 4 units, because 
2 


tim Ay = 2{[ tim — tim rim 
=2{ sas fash 


Since ultimate average problems are quite easy to compose, an instructor 
should have no trouble in giving his class a sufficient number of them. 


and 


| 
| 
A, 
1 n 
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ELEMENTARY PROBLEMS AND SOLUTIONS 
EpITED By Howarp EvEs, University of Maine 


Send all communications concerning Elementary Problems and Solutions to Howard 
Eves, Mathematics Department, University of Maine, Orono, Maine. This department 
welcomes problems believed to be new, and demanding no tools beyond those ordinarily fur- 
nished in the first two years of college mathematics. To facilitate their consideration, solutions 
should be submitted on separate, signed sheets, within three months after publication of prob- 
lems. 


PROBLEMS FOR SOLUTION 
E 1216. Proposed by N. A. Court, University of Oklahoma © 


The area of the triangle formed by the midpoints of three (not necessarily 
concurrent) cevians drawn through the three vertices of a given triangle is equal 
to one fourth of the area of the triangle determined by the feet of those cevians. 


E 1217. Proposed by Hiiseyin Demir, Zonguldak, Turkey 
Evaluate 


[J , 


d\n 
E 1218. Proposed by Robert Spira, Berkeley, Calif. 


Consider the two propositions: 
I. If (a, 6) =1, then ax+ 5 assumes infinitely many prime values. 
II. If (a, b) =1, then ax+5 assumes at least one prime value. 
I is Dirichlet’s theorem. Clearly I implies II. Show that II implies I. 


E 1219. Proposed by D. A. Robinson, University of Wisconsin 


The sequence {se}, where and Say =V. 2+-/s,, has a finite limit. 
Find this limit in closed form. 


E 1220. Proposed by Judith Blankfield, University of Illinois 
It is well known that if 
then 


> 


tan+1 


Fa 


Prove the following generalization which reduces to the above case when 6 =7. 


S r,/sin (6/2), 0 <@ < 2x. 


tmn+1 
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SOLUTIONS 
Total Number of Permutations of n Things 


E 1186 [1955, 655]. Proposed by David Freedman, Student, McGill University 


Show that the total number of permutations of m things is [n!e], where [x] 
denotes the greatest integer in x. 


Solution by P. B. Johnson, Occidental College. The total number P of permu- 
tations of m things is the sum of the number of permutations taken n, n—1, 
n—2,---,1,0 ata time. Hence 


P=) nl/r! = nle— >> 
r=0 


But 


0< n\/rl < 1)* = 1/n 1, 


ren+1 s=1 


whence the result. 

Also solved by A. N. Aheart, Ferrel Atkins, S. K. Berberian, D. G. Brennan, 
G. B. Charlesworth, John Christopher, A. E. Danese, R. E. Daniels and J. D. E. 
Konhauser (jointly), Marcel Delcourte, J. M. Elkin, Michael Goldberg, D. S. 
Greenstein, Nathaniel Grossman, M. J. Hellman, Vern Hoggatt, P. W. M. 
John, E. S. Keeping, D. C. B. Marsh, F. E. Milne, T. F. Mulcrone, C. S. 
Ogilvy, Hiram Paley, Walter Penney, C. F. Pinzka, L. A. Ringenberg, Azriel 
Rosenfeld, C. M. Sandwick, Sr., J. P. Scholz, W. A. Small, Bernard Smilowitz, 
O. E. Stanaitis, Art Steger, Chih-yi Wang, Alan Wayne, David Zeitlin, and the 


proposer. Late solutions by A. C. Goddard, A. R. Hyde, F. I. John, and M. S. 
Klamkin. 


A Limit Problem 


E 1187 [1955, 655]. Proposed by James Goode, Student, Georgia Institute of 
Technology 


Find 
2n 
lim >> 1/k. 
been 
I. Solution by G. B. Charlesworth, Hofstra College. By areas 
2n+1 2n 2n 
n kun n—1 
whence 


In (2m + 1)/n < y 1/k < In 2n/(n — 1). 
ben 


Zz 
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Therefore 


tn 
lim >> 1/k = In2. 


ken 


II. Solution by A. E. Danese, Mount Morris, N. Y. We have 
2n n—1 
> 1/k = 1/k 
ken kml k=l 


an) — (2 Inn) + 


Since the limit of each parenthesis is Euler’s constant y, we have 


2n 
lim >> 1/k = In 2. 


ken 


III. Solution by O. E. Stanaitis, St. Olaf College. From the definition of a 
definite integral 


tim = lim (1/n) $21/(1 + = f /(1 + 2) = In 2. 


IV. Solution by Albert Wilansky, Lehigh University. One may easily show by 
induction that 


k=l 


whence the required limit is In 2. 

Also solved by D. S. Adorno, W. A. Al-Salam, J. E. Darraugh, Marcel 
Delcourte, Hiiseyin Demir, I. A. Dodes, E. S. Eby, J. M. Elkin, Calvin Fore- 
man, A. M. Glicksman, Michael Goldberg, D. S. Greenstein, R. E. Greenwood, 
Cornelius Groenewoud, Nathaniel Grossman, M. J. Hellman, A. S. Hendler, 
T. W. Hildebrandt, Vern Hoggatt, W. H. Hooker, F. A. Homann, P. W. M. 
John, P. B. Johnson, L. M. Jordan, E. S. Keeping, M. A. Kirchberg, J. D. E. 
Konhauser, H. D. Lipsich, Kenneth Loewen, D. C. B. Marsh, Eugene McLach- 
lan, Elliott Mendelson, J. D. Miller, T. A. Newton, C. S. Ogilvy, Ingram Olkin, 
Hiram Paley, F. D. Parker, Walter Penney, W. J. Pervin, L. A. Ringenberg, 
Azriel Rosenfeld, J. W. Ross, C. M. Sandwick, Sr., L. P. Scholz, D. D. Strebe, 
Chih-yi Wang, R. M. Warten, Alan Wayne, David Zeitlin, and the proposer. 
Late solutions by A. R. Hyde, F. I. John, M. S. Klamkin, and C. F. Pinzka. 


An Orthocentric Tetrahedron 
E 1189 [1955, 655]. Proposed by N. A. Court, University of Oklahoma 


If two pairs of spheres with non-coplanar centers are such that each sphere 


| 

| 
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of one pair is orthogonal to the two spheres of the other pair, then the tetrahe- 
dron formed by the centers of similitude of the two pairs of spheres is orthocentric. 


Solution by P. D. Thomas, Air Force Armament Center, Eglin Air Force Base, 
Florida 


A consequence of the orthogonality conditions stated is that the two coaxial 
pencils determined by the two pairs of spheres as given are orthogonal. That is, 
every sphere of either pencil is orthogonal to every sphere of the other, and of 
course the two lines of centers are orthogonal. 

Now the sphere of similitude of each given pair (the sphere having for 
diameter the segment determined by the centers of similitude of the pair) be- 
longs to the coaxial pencil determined by that pair, and the two respective 
spheres of similitude are orthogonal. 

But the ends of two skew rectangular diameters of two orthogonal spheres 
are vertices of an orthocentric tetrahedron. Hence the proposition. 

See N. A. Court, Modern Pure Solid Geometry, pp. 153, 168, 179, 185, 189, 
269. 


Also solved by Roscoe Woods and the proposer. 
A Lucky Larry 
E 1190 [1955, 655]. Proposed by R. A. Miller, University of Mississippi 


A student wrote the following in solving the problem, “Find the value of 
(3+x)?/* when x =3/8:” 


(3 + x)?/* = (3 + 3/8)?/* = 2+ 1/4 = 9/4. 
Investigate the conditions under which 
(a + b/c)™!/* = (a + b/c)(m/n), 
where a, b, c, m, nm are positive integers. 


Solution by E. P. Starke, Rutgers University. First we note that two different 
positive integral k (k>1) powers must differ by more than &, since 
— s* = (¢ — s)(r>* + + + 


and r—s21 and the second expression in parentheses exceeds k. 
Consider 


= u(m/n), 


where u is positive and rational, and m, n are relatively prime positive integers. 
If m>n, we have u™-*=(m/n)*. Since m—n and n are relatively prime, both 
m and m must be perfect (m—n)th powers. But they differ by m—n. By our 
earlier remark this is impossible unless m—n =1. Hence m=n-+1 and 


w=a+ b/c = (1+ 1/n)* 
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with a arbitrary. 
The case m =n is trivial. If m<n, we have (n/m)"=u"-™ and, by the same 
argument, »=m+1 and 


u=a+t+ b/c = (1+ 1/m)™., 


The case cited by the Proposer results from n=3, m=2. 

Also solved by G. B. Charlesworth, Marcel Delcourte, I. A. Dodes, Michael 
Goldberg, Cornelius Groenewoud, Virginia Hanly, D. C. B. Marsh, O. C. Nuss- 
mann, C. S. Ogilvy, Walter Penney, Susan Pyeatt, Azriel Rosenfeld, G. J. 
Simmons, O. E. Stanaitis, Chih-yi Wang, R. M. Warten, H. R. Wesson, and 
Hazel Schoonmaker Wilson. Late solutions by A. R. Hyde, F. I. John, M. S. 
Klamkin, and M. A. Rashid. 


ADVANCED PROBLEMS AND SOLUTIONS 
EpitTep By E. P. STaRKE, Rutgers University 


Send all communications concerning Advanced Problems and Solutions to E. P. Starke, 
Rutgers University, New Brunswick, New Jersey. All manuscripts should be typewritten 
with double spacing and margins at least one inch wide. Problems containing results believed to 
be new or extensions of old results are especially sought. Proposers of problems should also en- 
close any solutions or information that will assist the editor. In general, problems in well 
known text books or results in readily accesstble sources should not be proposed for this de- 
partment. 


PROBLEMS FOR SOLUTION 
4688. Proposed by A. H. Clifford, Tulane University 


What positive integers n have the property that every positive integer less 
than m is expressible as the sum of distinct divisors of m? 


4689. Proposed by D. J. Newman, Republic Aviation Corporation, Farming- 
dale, N. Y. 


Let f(x) be any function such that f’’’(x) 20, f(n)~f(n+1). Prove that 


tends to as t—1-. 


4690. Proposed by Victor Thébault, Tennie, Sarthe, France 


Being given a tetrahedron ABCD and the tetrahedron A,B,C,D, obtained 
by passing planes through A, B, C, D parallel to the opposite faces of ABCD, 
show that 


PA* + PB* + PC? — 2PD* — PD} 


1 
i 
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is a constant independent of the position of point P. Extend this property to a 
skew polygon of vertices. 

4691. Proposed by R. E. Johnson, Smith College 


A weakly prime ideal of a ring R is any ideal J having the property that 
either aRCJ or RaCI implies that a is in J. Give an elementwise characteriza- 
tion of the unique minimal weakly prime ideal of R. 


4692. Proposed by A. E. Currier, United States Naval Academy 
Given f(x) defined by 


1 


prove that 
3+ 8/6 1 
1 — —— < liminf <-— 
(1) 750 im inf f(z) < 
1 3+ 8/6 
2 — <ili < 
(2) 750 ~ 
SOLUTIONS 
Digits of the Fibonacci Numbers 


4609 [1954, 646]. Proposed by Vern Hoggatt, San Jose State College, Cali- 
fornia 


Let gx(F,) be the kth digit (from the right) of F,, the nth Fibonacci number 
(Fo=0, Fi=1, Fase2=Faiit Fn). Let s be defined such that q(F,)=0 for 
n=0, 1, 2,-+-+,s but ge(F.4:) 0. (1) Show that the sequence g:(F,) has at 
least another run of s+1 consecutive zeros. (2) Prove or disprove: the sequence 
qx(F,) has never more than s+1 consecutive zeros. 


Editorial Note. No communication regarding this problem has been sent in, 
but some closely related questions are answered in the Proposer’s paper, A type 
of periodicity for Fibonacci Numbers, Mathematics Magazine, vol. 28, 1955, pp. 
139-142. 


Monic Irreducible Polynomials 
4624 [1955, 126]. Proposed by H. S. Shapiro, New York University 
Let J, denote the field of integers modulo p, where p is a prime, and J,[x] 


the ring of polynomials with coefficients in J,. Then the number of monic ir- 
reducible polynomials of degree n in J,[x] is equal to 


1 
— 
n 


din 


| 

= 
¢ 
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where yp is the Mobius function. 


Solution by B. M. Stewart, Michigan State University. A polynomial of degree 
s which is irreducible in J,(x) divides x*"—x if and only if s divides n. (Ex. 1, 
p. 180, MacDuffee, Introduction to Abstract Algebra.) Furthermore, x?"—x has 
no repeated factors. Hence if N(p, s) indicates the number of distinct monic 
irreducible polynomials of degree s in J,(x), then 


= sN(p, 5) 


summed over the positive divisors s of n. 
This is a perfect situation for applying the inversion formula which involves 
the Mdbius function yp. The result is 


nN(p, n) = 


summed over the positive divisors d of n. 

Also solved by P. T. Bateman, Robert Breusch, L. Carlitz, L. E. Clarke, 
Harley Flanders, S. W. Golomb, D. C. B. Marsh, Arnold Walfisz, Koichi 
Yamamoto, and the Proposer. 


Editorial Note. The result may be found in many places including Dedekind, 
Journ. reine u. angew. Math., vol. 54 (1857), pp. 19-22; L. E. Dickson, Linear 
Groups (1901), p. 18; H. Hasse, Zahlentheorie (1949) p. 32—references supplied 
by readers. 

Legendre Polynomials 

4628 [1955, 186]. Proposed by L. Carlitz, Duke University 

Let p be an odd prime and let P,,(x) denote the Legendre polynomial of 
degree m, where p=2m-+1. Show that 

0 od odd 
P.(3) = Pa(—3) = { (mod (m 
2a (mod (m even), 


where a is the unique odd integer determined by 
a@=b+1 (mod 4). 


Solution by the Proposer. Since P,,(—x)=(—1)"P,,(x), it suffices to discuss 
P,,(3). 

We shall make use of the following well known formula for the Legendre 
polynomial: 


(1) — 2x) = F(—m, m + 1; 1; 2), 


where F(a, b; c; x) denotes the hypergeometric function. (For proof, see Whit- 
taker and Watson, Modern Analysis, 4th edition, p. 311). 
Now it is clear from (1) that 


1 


| 
| 
| 


1956] ADVANCED PROBLEMS AND SOLUTIONS 349 


P,(1 — 2x) = (mod 


For x = —1 this becomes 
m 2 

(2) P.(3)= (mod 9). 
ka 


For m odd the sum in the right member of (2) evidently vanishes since the kth 
term and the (m—k)th term cancel. Hence P,,(3) =0 (mod ) in this case. 
When m is even, let m=2n; then we have 


3a 2n\? 2n 
(3) > = 

b=0 k n 
This may be proved by comparing coefficients of x in both members of the 
identity 


(x — + 1)* = (x? — 1)™. 
It follows from (2) and (3) that 


P»,(3) = (mod = + 1). 


But by a formula of Gauss (Theorie residuorum biquadraticorum, Werke, vol. 2, 
p. 90) 


= 2a (mod 9), 


where a is the unique odd integer determined by p=a?+5?, a=b+1 (mod 4). 
Since = (2/p) where (2/p) is the Legendre symbol, it 
follows that P:,(3)=2a (mod 


The following additional remark has some interest. If p =2m-+-1=3 (mod 4), 
then 


The problem was suggested by its connection with elliptic functions and 
particularly with the lemniscate function (see A. Hurwitz, Mathematische 
Annalen, vol. 51 (1899), pp. 196-226; Mathematische Werke, Basel, 1933, vol. 2, 
pp. 342-373). The following general result may be stated. Assume that the 
elliptic function sn x admits of complex multiplication so that the period quo- 
tient belongs to an imaginary quadratic field of discriminant d. Let k* denote 
the corresponding singular modulus and let » be an odd prime such that the 
Legendre symbol (d/p) = —1. Then P,,(1—2k*) =0 (mod p). See Mathematische 
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Annalen, vol. 127 (1954), pp. 162-169 and Monatshefte fiir Mathematik, vol. 58 
(1954), pp. 77-90. This result does not, however, determine the residue of 
P»(1—2k?) (mod p) when (d/p) = +1. 


Subring of a Radical Ring 


4629 [1955, 186]. Proposed by M. Henriksen and S. Perlis, Purdue Univer- 
sity 


A commutative ring R is called a radical ring (in the sense of Jacobsen) if 
every element has a quasi-inverse; 7.e., for every a in R, there is an x in R such 
that a+x+ax=0. In particular, if every element of R is nilpotent, then R is a 
radical ring. Need every subring of a radical ring be a radical ring? 

I. Solution by A. J. Goldman, Princeton University. We show first that an 
arbitrary commutative ring R (radical or not) without unity can be imbedded 
in a radical ring R’’. Imbed R in a commutative ring R’ with unity (it is well 
known that this can be done) and let F be the quotient field of R’. Then R is 
imbedded in the subring R” of F consisting of all elements of the form a/(1+5), 
where a and b range through R; it is clear that R” is actually a ring. R” is radical 
because a/(1+0) has the quasi-inverse —a/(1+a+5). 

To answer the question in the negative, it now suffices to exhibit a commuta- 
tive ring R without unity which is not radical. The even integers form such a 
ring, since the element 2 has no quasi-inverse. 

II. Solution by P. M. Cohn, the University, Manchester, England. Let R be 
the ring of formal power series in x over a field F, with zero constant term. Then 
R is a radical ring since any f in R has the quasi-inverse —f+f?—f*+ ---.The 
polynomials in R form a subring S which is not a radical ring since it contains 
x but no quasi-inverse of x. For if 


f = ax" +--+ + (a, 0) 


were a polynomial satisfying x+f+xf=0, then, by equating the coefficient of 
x"*+! to zero we obtain a; =0, which is a contradiction. Thus no non-zero element 
of S is a quasi-inverse of x; also zero is not. 

Also solved by W. R. Abel and C. D. Gorman, E. H. Batho, S. W. Golomb, 
C. R. Phelps, Gustave Rabson, R. B. Reisel, D. W. Sasser, W. R. Scott, M. F. 
Smiley, M. L. Tomber, Jerrold Yos, and the Proposer. 


A Continued Fraction Related to the Bessel’s Functions 


4631 [1955, 187]. Proposed by F. H. Northover, Memorial University of 
Newfoundland 


Evaluate the unending continued fraction 


Solution by D. J. Newman, Republic Aviation Corporation, Farmingdale, 


__] 
| 
| 
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N. Y. Applying the method of continued fractions to the differential equation 
y=xy"’ (See Ince, Differential Equations, p. 178) we obtain 


for some solution y. (Namely, letting x0, it is that solution for which y(0) =0, 
y’(0) =1). All that remains then is to solve the equation 


y= xy”, 0)=0, y(0)=1 


and our answer will be y’(1)/y(1). 

The differential equation is of the Bessel type and the required solution is 
then y=x"2J,(2x/*), whence y’=J(2x"/*). This gives as our final answer 
Io(2)/Ih(2), or in series form 


2 4 1 
nin + 1)! 


A. numerical value, 1.433127, is easy to obtain either from the series or from 
convergents to the continued fraction. 

Also solved by Leonard Carlitz, M. Delcourte, V. C. Harris, P. G. Kirmser, 
D. H. Lehmer, Erich Michalup, C. D. Olds, C. F. Pinzka, C. C. Sams, R. E. 
Shafer, M. R. Spiegel, Chih-yi Wang, Alan Wayne, and the Proposer. 


Editorial Note. The proposed continued fraction is a special case of many in 
the literature. Solvers referred to Chrystal, Algebra, vol. 2, p. 494, and pp. 520— 
522; Erdelyi, et al., Higher Transcendental Functions, vol. 2, formula 7.11(23); 
H.S. Wall, Continued Fractions, p. 349; Perron, Die Lehre von der Kettenbriichen, 
p. 299 and pp. 479-480. 


RECENT PUBLICATIONS 
Epitep By E. P. VANcE, Oberlin College 
All books for review should be sent directly to E. P. Vance, Oberlin College, Oberlin, Ohio. 


Transform Calculus with an Introduction to Complex Variables by E. J. Scott. 
Harper and Brothers, New York, 1955. 8+330 pages. $7.50. 


In the introduction the author indicates that this book is directed toward 
advanced undergraduate and graduate students, mainly in physics, chemistry, 
and engineering. For this audience it is excellent, providing them with many 
techniques and specific results of much importance in their work. Mathematical 
rigor is attempted within reason, but the author’s emphasis is more on vigor 
than on rigor. He does not attempt the most general known statement of a 
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theorem if its proof would take him through circuitous pedantics. Thus he con- 
fines himself to piecewise continuity in his statement of the Fourier integral 
theorem (which he does not prove). Some theorems are stated imperfectly, and 
the English is not always perfect. Some discussions might be written more satis- 
factorily from the logician’s viewpoint, but it would be difficult to do this with- 
out spoiling their pedagogic value. There are some misprints. 

Chapter I gives the highlights of the theory of functions of a complex vari- 
able. Chapter II introduces the Fourier and Laplace transforms and their in- 
verses. Some highlights of the theory of the Laplace transform are given in 
Chapter III, and are applied to ordinary differential equations with constant 
coefficients in Chapter IV. Chapter IV contains a number of examples of prac- 
tical interest from electric circuit theory and elasticity. 

In Chapter V matrix concepts are introduced and applied to systems of 
linear differential equations of the type that occur in multiple degree of freedom 
electrical and mechanical systems. These are solved using the Laplace trans- 
form. A number of examples involving actual applications to electrical and 
mechanical networks are given. 

In Chapter VI, the Laplace transform is applied to the solution of difference 
equations; in particular séveral applications are made to mechanical and electri- 
cal problems. 

Chapter VII is a long chapter in which the Laplace transform is applied to 
the solution of linear partial differential equations. Diffusion is treated in detail, 
and some transmission line problems are solved. Longitudinal and transverse 
vibrations of a bar are treated, as well as transverse vibrations in a string. The 
emphasis is upon the particular applications. Nothing like a unified theory of 
elliptic, hyperbolic, or parabolic equations is attempted. 

In Chapter VIII, Gamma, Error, Fresnel, and Bessel functions are briefly 
introduced, as well as asymptotic series and some simple integral equations. 
The author is to be commended for employing the standard notation Y,(¢) for 
the now universally employed second solution of the Bessel equation which is 
due to Weber. Jahnke and Emde erroneously attributed this function to Neu- 
mann and changed the notation to N,(¢). Unfortunately many electrical engi- 
neers in this country have followed this schismatic practice. 

Chapter IX treats some partial differential equations of a more complicated 
nature which involve the special functions introduced in Chapter VIII. A 
number of problems of practical interest are worked. 

In the last chapter, Chapter X, the author introduces what he calls finite 
Fourier transforms. The resulting theory is equivalent to the theory of Fourier 
series. He also introduces what he calls the finite Hankel transform. The result- 
ing theory is equivalent to the theory of Fourier-Bessel series. 

This book will provide engineers and physical scientists with a good many 
techniques that they will find important. It will not be satisfactory as a last 
word on the subjects treated to the mathematics major (for instance, a mathe- 
matics major ought to be able to prove Fourier’s theorem). On the other hand, 
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a good many mathematics students would be much better off for having had a 
course of this type. There might then be fewer members of that forlorn group 
of non-productive Ph.D.’s who know theorems precisely, but have no idea what 
in the world to de. with them. 
A valuable feature of the book is its many examples and illustrations. The 
typography is excellent. 
E. PINNEY 
University of California 
Berkeley, California 


Numerical Methods. By A. D. Booth. New York, Academic Press. 1955. 7+195 
pages. $6.00. 


This is a small book (187 pages) covering a wide range of topics in Numerical 
Analysis. Since it covers such a wide range of topics the treatment of each sub- 
ject is necessarily brief. The author states that the book was prepared with 
automatic computing equipment in mind; therefore, there is emphasis on meth- 
ods suitable for use with automatic computing equipment. 

The book starts with interpolation formulas using finite differences. Numer- 
ical integration and differentiation is the subject of the next chapter. Summation 
of series and the solution of ordinary and partial differential equations, as well 
as the solution of linear and non-linear systems of equations, are treated. The 
book ends with approximating functions, Fourier synthesis and analysis, and 
Monte Carlo methods. 

There are some misprints in the book. For example, on Page 2 read SWAC 
for SEAC. Clearly, formula 4.6.12 on Page 39 should read s,y"/m instead of 
Say”. In many places in Chapter 6 the error terms of formulas involve factors 
such as f‘ max. These should be replaced by terms of the form f(s) where s 
is some point in the interval under consideration. With these indications of lack 
of care in the statement of formulas, one realizes that other formulas in the book 
must be used with caution. 

However, the development of high speed computing equipment has increased 
the demand for material on numerical analysis so this book should be part of the 
library of all those interested in computation. 

H. D. HusKEy 
University of California 


Automatic Feedback Control System Synthesis. By J. G. Truxal. McGraw-Hill, 
New York. Electrical and Electronic Engineering Series, 1955, xiii+-675 
pages. $12.50. 


This book is intended as an engineering text and reference. It achieves the 
second aim much more successfully than the first. The book is based on a second 
course in servomechanisms as taught by the author. The book contains many 
cases where initial references are made apparently to content material of the 
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first course. In some of these instances a text reference is given, but this pro- 
cedure is confusing when reading the book per se. 

The book is replete with completely solved examples. In fact there is a tend- 
ency to explain by using many specific illustrations. The book might have been 
considerably reduced in size with little loss if proofs were given or generalized 
illustrations used. The detail with which examples are worked out is also ex- 
cessive. At this level for example, it should not be necessary to work out poly- 
nomial long divisions in detail. 

The contents are divided as follows: 

Chapter 1, Introduction; develops Laplace Transform Techniques. 

Chapter 2, Signal-flow Diagrams and Feedback Theory; and Chapter 3, 
Synthesis of RC Networks; develop the basic theory of network synthesis 
and circuit stability. 

Chapter 4, Root-locus Methods; Chapter 5, Synthesis through Pole-zero 
Configurations; and Chapter 6, Design in the s Plane; constitute the major 
contribution of the book. Here the details of techniques for linear network 
synthesis are described fully. 

Chapter 7, Principles of Statistical Design; Chapter 8, Application of 
Statistical Design Principles; and Chapter 9, Sample-data Control Systems; 
develop the statistical treatment leading to system design techniques where 
noise is involved. — 

Chapter 10, Non-linear Systems and Describing-function Analysis; and 
Chapter 11, Phase-plane Analysis; discuss several methods of treating non- 
linear systems. 

From a mathematical point of view it is distressing to find the author ac- 
centuating Laplace Transform theory as the only essential mathematical tool 
necessary for servo-engineers. A much healthier approach would be to accentu- 
ate the value of a strong mathematical background for all engineers. 

Generally the book is well written and well printed. Few misprints were 
noted. The book is comprehensive and should prove useful as a reference for the 
field. 

Brock 
Electro Data Corporation 
Pasadena, California 


NEW BOOKS RECEIVED 


Modern Trigonometry. By W. A. Turledge and J. A. Pond. Englewood Cliffs, 
N., J., 1956. 11+243 pages. 

An Introduction to Mathematics. By L. E. Boyer. New York, Henry Holt 
and Co., 1956. 16+528 pages. $5.25. 

Nellie Landblom's Copybook for Beginners in Research Work. By Nellie 
Thompson Landblom. Multigraphed book available from author, P. O. Box 
248, Fort Collins, Colorado, 1955. 118 pages. $2.95. 
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Elements of Business Mathematics for Colleges. Llewellyn R. Snyder. New 
York, McGraw-Hill Book Co., 1956. x +249 pages. $3.75. 

Carl Friedrich Gauss: Titan of Science. A Study of His Life and Work. By 
G. W. Dunnington. New York, Exposition Press, 386 Fourth Ave., New York 
16, N. Y., 1955. xi+479 pages. $6.00. 

Surface Area. By Lamberto Cesari. Princeton, N. J., Princeton University 
Press, 1956. x +595 pages. $8.50. 

Articles et Memoires—I, La variable complexe, II, Le champ reel—Notices— 
Arnaud Denjoy. Paris, Gauthier-Villars, 1955. Vol. 1, x +507 pages and Vol. 2, 
vi+600 pages. 2 volumes $15.03. 

Calculus, Differential and Integral, with Problems and Solutions. By G. M. 
Peterson and R. F. Graesser. Ames, Iowa. Littlefield, Adams and Co., 1956. 
10+321 pages. $1.75. 

Through the Mathescope. By C. S. Ogilvy. New York, Oxford University 
Press. 1956. viii+162 pages. $4.00. 

Elementary Topology. By D. W. Hall and G. L. Spencer. New York, John 
Wiley and Sons, 1955. xii+303 pages. 

Tracts for Computers, No. XX VI, Correlated Random Normal Deviates. By 
E. C. Fieller, T. Lewis and E. S. Pearson. New York, Cambridge University 
Press, 1956. 60 pages. $2.00. 

Spheroidal Wave Functions, Including Tables of Separation Constants and 
Coefficients. By J. A. Stratton, P. M. Morse, L. J. Chu, J. D. C. Little and F. J. 
Corbato. Published jointly by the Technology Press of M.I.T. and John Wiley 
and Sons, Inc., New York, 1956. xiii+-613 pages, $12.50. 

Plane Waves and Spherical Means Applied to Partial Differential Equations. 
By Fritz John. New York, Interscience Publishers, Inc., 1955. 180 pages. $4.50. 

Experimental Design. By W. T. Federer. New York, The Macmillan Com- 
pany, 1955. xix+544+-47 pages of problems. $11.00. 

Mathematics of Finance, Second Edition. By P. M. Hummel and C. L. See- 
beck, Jr. New York, McGraw-Hill Book Company, 1956. xi: 372 pages. $4.75. 

Analytic Geometry, Fifth Edition. By C. E. Love and E. D. Rainville. New 
York, The Macmillan Co., 1955. xiv+302 pages. $4.00. 

Elementary Statistics for Students of Social Science and Business. By R. C. 
Sprowls. New York, McGraw-Hill Book Company, 1955. xiii+392 pages. $5.50. 

Analytic Geometry and Calculus. By T. S. Peterson. New York, Harper and 
Brothers, 1956. ix+456 pages. $5.50. 

Introduction to Numerical Analysis (International Series in Pure and Applied 
Mathematics). By F. B. Hildebrand. New York, McGraw-Hill Book Company, 
1956. x +511 pages. $8.50. 

Mathematical Theory of Elasticity, Second Edition. By I. S. Sokolnikoff. 
New York, McGraw-Hill Book Company, 1956. xi+476 pages. $9.50. 

Plane Trigonometry, Third Edition. By A. L. Nelson and K. W. Folley. 
New York, Harper and Brothers, 1956. xi+134 pages. $3.50. 
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Numerical Analysis. By Zdenek Kopal. New York, John Wiley and Sons, 
Inc., 1955. xiv-+556 pages. $12.00. 

Methods in Numerical Analysis. By K. L. Nielsen. New York, The Macmil- 
lan Company, 1956. viii+382 pages. $6.90. 

An Introduction to Linear Algebra. By L. Mirsky. New York, Oxford Uni- 
versity Press, 1956. xii+433 pages. $5.60. 


NEWS AND NOTICES 
EpiTep By Epita R. SCHNECKENBURGER, University of Buffalo 


Readers are invited to contribute to the general interest of this department by sending news 
items to Edith R. Schneckenburger, University of Buffalo, Buffalo 14, New York. Items 
should be submitted at least two months before publication can take place. 


SUMMER SESSIONS 


The following institutions announce advanced courses in mathematics for 
the summer of 1956. 

Michigan State University. June 26 to August 3: Professor Arnold, correlation 
analysis; Professor Hill, advanced calculus I, functions of a complex variable; 
Professor J. B. Kelly, differential equations, matrices and groups; Professor 
Meixner, potential theory, mathematical theory of continuous matter; Professor 
Stewart, projective geometry. June 26 to August 24: Dr. Berberian, mathe- 
matics of engineering; Professor Blair, functions of a real variable, algebraic 
systems of continuous functions; Dr. Faith, advanced calculus II; Professor 
Hannan, analysis of variance; Professor L. Kelly, differential equations, vector 
analysis, combinatorial topology; Mr. Laurent, statistical methods in engineer- 
ing. 
Stanford University. June 25 to August 17: Professor Weinberger (University 
of Maryland), approximation of quadratic functionals; Professor Lax (New 
York University), selected topics from analysis; Professor Spencer (Princeton 
University), seminar. 

State University of Iowa. June 12 to August 8: Professor Hogg, introduction 
to mathematical statistics, statistical hypotheses; Professor Johnson, astron- 
omy; Professor Knowler, operations research, actuarial aspects of casualty insur- 
ance; Professor Muhly, advanced calculus, group theory; Professor Price, the 
supervision of mathematics, laboratory in the construction of teaching aids in 
mathematics; Professor Smiley, matrices and determinants, introduction to 
advanced mathematics; Professor Woods, studies in secondary geometry, differ- 
ential equations. 
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University of Michigan, School of Education. June 25 to July 13: The Work- 
shop for College Professors will present a program designed to help college teach- 
ers develop new ideas and fresh materials for their academic courses. Such topics 
as course planning, teaching techniques, and evaluation will be discussed. For 
further information write to Director A. D. Henderson, 2442 U.E.S., University 
of Michigan, Ann Arbor, Michigan. 

University of Minnesota, College of Science, Literature, and the Arts. June 11 
to July 14: Professor Carlson, differential equations, non-euclidean geometry; 
Dr. Orey, Fourier series and orthogonal functions; Visiting Professor Roy, 
analysis of variance and multivariate analysis II, statistical models in the social 
sciences. July 16 to August 18: Professor Gelbaum, integral equations, inter- 
mediate differential equations; Professor Hatfield, Jr., Fourier series and orthog- 
onal functions, II. June 18 to July 21: A Summer Institute for High School 
Teachers of Mathematics, carrying graduate credit, will be given with the sup- 
port of the Hill Family Foundation. 

University of Pittsburgh. June 11 to July 20, July 23 to August 31: Professors 
Blumberg and Christiano, differential equations, advanced calculus; Professor 
Bryson, vector analysis; Professor Taylor, functions of a complex variable; 
Professor Levine, modern algebraic theories; Professor Laush, infinite series, 
Lie theory of groups (second six weeks); Professor Bompiani, geometry of 
differential elements; Professor Elyash, advanced function theory (first six 
weeks). June 18 to August 10 (evenings): Professor Knipp, mathematics of 
modern engineering; Professor Bryson, Laplace transform theory and applica- 
tions (second half of two semester course). July 2 to August 10: Professor Vog- 
eley, differential equations for students of engineering; Professor Edwards, 
mathematical theory of statistics; Professor Teats, history of mathematics; Pro- 
fessor Myers, recreational mathematics for teachers, theory of equations; Mr. 
Kachun, navigation for teachers; Professor Knipp, teaching of secondary school 
mathematics. 

University of Virginia. June 18 to August 14: Professor Ball, advanced cal- 
culus and applied mathematics; calculus of variations; Professor Floyd, topo- 
logical groups. July 3 to August 14: Professor Ball, foundations of geometry. 


CLIFFORD BREWSTER UPTON FELLOWSHIP 


The American Book Company has established a fellowship in honor of Pro- 
fessor Emeritus Clifford Brewster Upton for graduate study in the Teaching of 
Mathematics at Teachers College, Columbia University. Three thousand dol- 
lars will be awarded annually to a student in the teaching of mathematics in 
elementary schools, in secondary schools, or in teacher training institutions. 
Inquiries should be made to the Chairman of the Committee on Fellowships 
and Scholarships, Teachers College, Columbia University, New York 27, New 
York, before May 1, 1956. Formal applications, together with all supporting 
materials, must be filed with the Committee by June 1, 1956. 
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CALIFORNIA CONFERENCE FOR TEACHERS OF MATHEMATICS 


The 1956 California Conference for Teachers of Mathematics will be held 
at the University of California at Los Angeles from June 20 through July 3. The 
purpose of this conference is to bring together teachers interested in mathemat- 
ics—arithmetic through calculus—to study common problems in the teaching 
of mathematics. For further information write to University of California Ex- 
tension, Mathematics, Los Angeles 24, California. 


PERSONAL ITEMS 


Assistant Professor J. W. Lindsay of Texas Technological College has been 
awarded a Danforth Foundation Teachers Study Grant for the year 1956-1957. 
He will spend the year at the University of California at Los Angeles. 

Illinois Institute of Technology announces the following: Mr. Michel Mc- 
Kiernan, previously a graduate assistant at the Institute, and Dr. Abe Sklar, 
formerly an assistant at California Institute of Technology, have been appointed 
to instructorships. 

Kansas State College reports: Professor P. M. Young has been appointed 
Associate Dean of the School of Arts and Sciences; Professor R. G. Sanger, 
formerly the Acting Associate Dean of the School of Arts and Sciences, has re- 
sumed his duties as Head of the Department of Mathematics. 

At Wartburg College: Professor J. O. Chellevold, head of the Department of 
Mathematics, has been appointed Acting Dean of the Faculty; Mr. Richard 
Wideranders of Estherville Junior College, Iowa, has been appointed to an 
instructorship. 

Mr. R. G. Albert, formerly a graduate student at Harvard University, has a 
position as a programmer for the Datamatic Corporation, Newton Highlands, 
Massachusetts. 

Mr. P. N. Armstrong, previously a technical engineer for the International 
Business Machines Corporation, Endicott, New York, is an engineer for the 
Magnavox Company, Santa Monica, California. 

Mr. R. V. Benson of Los Angeles City College has been appointed to an 
instructorship at Long Beach City College. 

Dr. W. J. Berger, formerly a mathematician at R.C.A. Service Company, 
Patrick Air Force Base, Cocoa, Florida, now has a position as a mathematician 
with the Raydist Navigation Corporation, Hampton, Virginia. 

Assistant Professor D. R. Bey of Illinois State Normal University has been 
promoted to an associate professorship. 

Dr. E. E. Blanche, previously an analyst with the Research and Develop- 
ment Division, Department of the Army, Washington, D. C., is President and 
Senior Research Scientist of Ernest E. Blanche and Associates, Rockville, 
Maryland. 

Mr. H. L. Carlson, formerly an electronic applications specialist for the 
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National Cash Register Company, Dayton, Ohio, is employed as a computing 
engineer by Douglas Aircraft Company, Santa Monica, California. 

Mr. D. A. Cope, recently an engineer for the Glenn L. Martin Company, 
Baltimore, Maryland, has accepted a position as a senior research engineer for 
the Consolidated-Vultee Aircraft Company, San Diego, California. 

Mr. H. N. Garber, formerly National Science Foundation fellow at Massa- 
chusetts Institute of Technology, has a position as a management assistant with 
the Radio Corporation of America, Camden, New Jersey. 

Miss Jane C. Ingersoll, previously an operations analyst at the Operations 
Research Office, Chevy Chase, Maryland, has been appointed to an instructor- 
ship at the University of Maryland. 

Mr. Seymour Landau, formerly an instructor at the University of Connecti- 
cut, is now an actuary for the New York City Transit Authority, New York. 

Mr. R. L. Mentzer of Teachers College of Connecticut has been promoted to 
an assistant professorship. 

Associate Professor E. N. Nilson of Trinity College, Connecticut, has a posi- 
tion as a project engineer for Pratt and Whitney Aircraft Corporation, East 
Hartford, Connecticut. 

Dean Emeritus and Professor J. R. Overman of Bowling Green State Uni- 
versity is retiring in June 1956. He was the first member of the teaching staff, 
first Chairman of the Department of Mathematics, and first Dean of the College 
of Liberal Arts, Dean of Men, and Dean of Faculties. The University’s Board of 
Trustees recently honored Dr. Overman by naming the new science building 
Overman Hall “in recognition of his unique service and contribution in the 
development and progress of the University.” 

Mr. J. L. Pallone, assistant actuary for Woodward and Fondiller, Inc., New 
York, New York, has been appointed to an instructorship at Polytechnic Insti- 
tute of Brooklyn. 

Assistant Professor H. C. Parrish of North Texas State College has been 
promoted to an associate professorship. 

Mr. J. B. Roberts of Reed College has been promoted to an assistant pro- 
fessorship. 

Dr. D. T. Sigley, formerly a physicist at Johns Hopkins University Physics 
Laboratory, is employed as Associate Director of the Engineering Laboratories, 
American Machine and Foundry Company, Greenwich, Connecticut. 

Assistant Professor James C. Smith, Jr., of Loyola University of Los Angeles 
is a teacher in the Extension Division, University of California at Los Angeles. 

Mr. G. E. Waldman, a teacher at Tottenville High School, Staten Island, 
New York, has been appointed Acting Chairman of the Department of Mathe- 
matics. 

Dr. E. A. Walker, recently a graduate student at the University of Kansas, 
is employed as a mathematician at the National Security Agency, Washington, 
D. C. 
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Dr. G. P. Weeg, a mathematician at Engineering Research Associates, St. 
Paul, Minnesota, is now employed in the Systems Analysis Department, Rem- 
ington Rand Univac Division of the Sperry Rand Corporation, St. Paul, Minne- 
sota. 

Mrs. Elizabeth S. Wolf, formerly an instructor at Indiana Technical College, 
has a position as a computer in the Advanced Development Division of the Avco 
Manufacturing Company, Stratford, Connecticut. 

Dr. Albert Wolinsky, previously an associate physicist with the Farrand 
Optical Company, New York, New York, has a position as a staff member of the 
General Precision Laboratory, Pleasantville, New York. 


Associate Professor Willie E. Clark, chairman of the Department of Mathe- 
matics of Agricultural, Mechanical, and Normal College, Pine Bluff, Arkansas, 
died on December 4, 1954. 

Captain C. E. Dimick, United States Coast Guard (ret.), professor and head 
of the Department of Mathematics of the United States Coast Guard Academy, 
died on January 2, 1956. He was a charter member of the Association. 

Professor Emeritus C. C. Grove of City College of the City of New York died 
on January 12, 1956. He was a charter member of the Association. 


THE MATHEMATICAL ASSOCIATION OF AMERICA 


Official Reports and Communications 
THE NOVEMBER MEETING OF THE PHILADELPHIA SECTION 


The annual meeting of the Philadelphia Section of the Mathematical Asso- 
ciation of America was held at the University of Pennsylvania on November 26, 
1955. Professor T. L. Koehler, Chairman of the Section, presided. There were 
87 present, including 64 members of the Association. By-Laws for the Section 
were adopted by vote of those present. 

Officers of the Section for the year 1955-1956 are: Chairman, Professor J. C. 
Oxtoby, Bryn Mawr College; Secretary, Professor G. C. Webber, University of 
Delaware; Members at Large—for three years, Professor N. J. Fine, University 
of Pennsylvania; two years, Professor T. L. Koehler, Muhlenberg College; one 
year, Dr. Bryant Tuckerman, Institute for Advanced Study. 

The following papers were presented: 


1. A report on the Ford Foundation study on the integration of high school and 
college mathematics, by Professor H. W. Brinkmann, Swarthmore College. 


Professor Brinkmann reported on the present status of the School and College Study of Ad- 
mission with Advanced Standing as it concerns mathematics. He gave a brief summary of the 
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results of the mathematics examinations given last spring and of the Conference of Mathematics 
Teachers held at Williams College in June. The work of the Committee on Mathematics, as well as 
that of the other committees, is being carried on by the College Entrance Examination Board. 
Examinations for advanced standing will therefore be available (beginning next May) for any 
candidates who wish to take advantage of them. 


2. Dedekind sums and classes of modular substitutions, by Professor Hans 
Rademacher, University of Pennsylvania. 


The modular matrices 
b 
1) @ b, integers, od — bc = 1). 


form a group I’. Two matrices M and M’ are called similar if there exists an LET such that 
M' = L“ML, 
The relation of similarity splits T into classes. An invariant for similarity is obviously the trace 
s=a+d. Another invariant can be defined as follows: 
b 
at forc = 0 


— 12sign -s(a, | ¢|) — 3 sign +4)) force #0 


= 


where s(h, k) is the Dedekind sum 
(hy hus 1 
The invariant W has the properties 


¥(M") = ¥(M). 


The inhomogeneous modular group I* (i.e. the factor group I'/E, where E consists of two elements, 
the unit matrix J and —J) can be generated by two elements T, U with 


0-1 i—1 
T= U = 
( 0) 
T? = = 1, 


Any MeEI* has therefore a unique representation 
M = + 


such that 


with ¢€0, é4:=0, +1, +1, 2, +--+». The paper gives then another determination of ¥(M) 
by means of the «. 


3. Flexagons, by Professor R. J. Wisner, Haverford College. 


During the last few years, a number of people have shown interest in certain geometric objects 
usually called flexagons. Three notes have been published about two of these: the one with three 
and the one with six physical faces. A method of constructing flexagons with » 23 physical faces, 
some properties of “regular” flexagons, a possible algebraic definition, and the use of flexagons as 
realizations of certain groups are shown. The results are a part of joint work by Professor C. O. 
Oakley and the speaker. 
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4. On differential operators, by Professor William Feller, Princeton Univer- 
sity. 


The differential operator Q defined by Qf =af’’+bf’+cf, where a>0, has the characteristic 
property: if f has a local minimum at xo and f(xo) =0 then Qf(xo) 20. The speaker described the 
form of the most general operator with this property, and the advantages arising from such 
generalized differential operators instead of Q. 


5. Pea soup, tripe and mathematics, by Professor Morris Kline, New York 
University. 


The speaker contends that the usual material taught to freshmen, namely, college algebra and 
trigonometry, is meaningless, unmotivated, lacking in esthetic qualities and incoherent. These 
deficiencies account for the usual negative student reaction to freshman mathematics. This paper 
also presents objections to the replacement of this conventional material by rigorous abstract 
mathematics such as set theory, symbolic logic, Dedekind cuts, groups and fields and the like. The 
recommendation made in this talk is that mathematics be tied intimately to physical problems and 
to its cultural setting and significance. Physical problems should provide the motivation, meaning 
and significance at the same time that the role of mathematics in idealizing and reasoning about 
physical problems is given the chief attention so that the course is definitely one in mathematics 
and not a part of a physical science course. Other values of mathematics that can be presented in 


such an approach are included. 


G. C. WEBBER, Secretary 


CALENDAR OF FUTURE MEETINGS 


Thirty-seventh Summer Meeting, University of Washington, Seattle, Wash- 


ington, August 20-21, 1956. 


Fortieth Annual Meeting, University of Rochester, Rochester, New York, 


December 29, 1956. 


The following is a list of the Sections of the Association with dates of future 
meetings so far as they have been reported to the Asseciate Secretary. 


ALLEGHENY MouNTAIN 


Eastern Illinois State College, 
Charleston, May 11-12, 1956. 

INDIANA, Wabash College, Crawfordsville, 
May 5, 1956. 

Iowa 

KANSAS 

KENTUCKY 


MARYLAND-DIstTRICT OF COLUMBIA-VIRGINIA, 
United States Naval Academy, Annapolis, 
Maryland, May 5, 1956. 

METROPOLITAN NEw YORK 


MICHIGAN 

Minnesota, Augsburg College, Minneapolis, 
May 5, 1956. 

MIssouRI 

NEBRASKA 


NORTHEASTERN, University of Connecticut, 
Storrs, November 24, 1956. 

NORTHERN CALIFORNIA 

OHIO 

OxtaHoMA, Oklahoma City University, Oc- 
tober 26, 1956. 

Paciric NortHwest, Oregon State College, 
Corvallis, June, 1957. 

PHILADELPHIA 

Rocky Mountain, University of Utah, Salt 
Lake City, May 4-5, 1956. 

SOUTHEASTERN 

SOUTHERN CALIFORNIA 

SOUTHWESTERN 

TEXAS 

NEw York STATE 

Wisconsin, Marquette University, Milwaukee, 
May 12, 1956. 


FOR THE MATHEMATICIAN 
who’s ahead of his time 


1BM is looking for a special kind of mathe- 
matician, and will pay especially well for 
his abilities. 

This man is a pioneer, an educator—with a 
major or graduate degree in Mathematics, 
Physics, or Engineering with Applied 
Mathematics equivalent. 


You may be the man. 


® If you can qualify, you'll work as a special 
BE2i representative of IBM’s Applied Science 
S Division, as a top-level consultant to busi- 
een ness executives, government officials and 


scientists. It is an exciting position, crammed 
with interest, and responsibility. 


xX Employment assignment can probably be 
made in almost any major U. S. city you 
choose. Excellent working conditions and 


= employee-benefit program. 
a For applicants with the same basic qualifi- 
+ cations, opportunities are available to teach 


in this exciting, new field. 


Your reply will, of course, be held in the 
Strictest confidence. Write, giving full de- 
tails of education and experience, to: 


Dr. C. R. DeCarlo 


DIRECTOR, APPLIED SCIENCE DIVISION 


INTERNATIONAL BUSINESS MACHINES CORP. 
590 Madison Avenue, New York 22, N.Y. 


Hundreds of IBM electronic data processing ma- 
chines are already in use, and many more will be 
installed during 1956. These data processing ma- 
chines have become essential for computations in 
science, engineering, and business management. 


Producer of electronic 

data processing machines, 
electric typewriters, 

and electronic time equipment. 


: i on 
>| 
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RAYMOND W. BRINK’S Standard Textbooks 
A FIRST YEAR OF COLLEGE MATHEMATICS, 2nd Edition 


This book presents a complete and integrated course in college algebra, 
trigonometry, and analytic geometry together with a review, in the Ap- 
pendix, of elementary algebra. Exercises are graded in order of difficulty. 

725 pages, $5.00 


COLLEGE ALGEBRA, 2nd Edition 


Notable for its accuracy, logic, and flexibility, this first-year text presents 
a complete course in college algebra introduced by a thorough and 
systematic review of high-school algebra. 495 pages, $4.00 


INTERMEDIATE ALGEBRA, 2nd Edition 


Presenting in full the first fifteen chapters of College Algebra, this text 
is designed for college students who have had only one year of high- 
school algebra and who desire to prepare for college mathematics or 
science. 295 pages, $3.25 


APPLETON-CENTURY-CROFTS, INC. 
35 West 32nd Street New York 1, N.Y. 


For Analysis Group of expanding research 
and development laboratory. For those 
who are interested in research dealing with 


M at he matic ian * weapons systems development, peacetime 


applications of atomic energy, automation 


and data processing, two excellent 
| openings are available: 


- PhD or equivalent, with training or experience in OPERATIONS RESEARCH. 
Familiarity with probability theory, linear programming, game theory, 
information theory, optimization procedures, and 
other OR techniques, very desirable. 


+ Masters degree. To investigate problems in ballistics, aerodynamics, optics and 
other fields of applied mathematics. May also conduet operational 
analyses on missile, fire control, navigational and other systems. Some 
familiarity with digital computer programming desirable. 


These openings require men with vision and initiative. Our modern laboratory 
provides a professional working atmosphere and the location in a 

quiet suburban area provides pleasant living and working with easy access to the 
cultural and educational facilities of New York City. 


All inquiries in confidence. Please send resume including salary desired, to 
MR. FRED SCHWARZ 


VITRO LABORATORIES vivo corn ot 


200 Pleasant Valley Way, West Orange, N.J. 


< 
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For summer and fall classos ... 


THEORY OF NUMBERS 
By B. M. Stewart, Michigan State University 


“,.. an interesting choice of subject matter. . .. The student using it would be 
left in no doubt as to the importance of the topic in hand, nor of the existence 
of subtleties in the proofs, This is a highly desirable characteristic in a text for 
a beginning theorist.”—Julia Wells, Connecticut College. 


1952 261 pages $5.75 


A SURVEY OF MODERN ALGEBRA, Revised Edition 


By Garrett Birkhoff, Harvard University, and Saunders MacLane, Uni- 
versity of Chicago 


“The manner in which familiar illustrations accompany new definitions and 
concepts is commendable. The authors have included a wide variety of exercises 
which should serve the desires of different students. The topic of matrices has 
been treated very well.” —J. Ray Hanna, University of Wichita 


1953 * 472 pages « $6.50 


A BRIEF SURVEY OF MODERN ALGEBRA 
By Garrett Birkhoff and Saunders MacLane 


The authors have converted the first ten chapters of the revised Survey of 
Modern Algebra into a brief survey which is suitable for shorter courses in linear 
algebra or modern algebra. It provides a clear and adequate introduction to the 
main ideas of modern algebra, omitting certain advanced topics. 


1953 ° 276 pages ° $4.75 


ELEMENTARY THEORY OF EQUATIONS 
By Samuel Borofsky, Brooklyn College 


“I have found my choice of texts in Elementary Theory of Equations . . 
The book has been composed with a style. The author seems to be speaking to 
the readers in a conversational tone. . . . The material is just about right for an 
introductory course to undergraduates.”—Norman Dodson, Lenoir Rhyne College. 


1950 ° 302 pages ° $4.75 


VECTOR AND TENSOR ANALYSIS 
By Nathaniel Coburn, University of Michigan 


This unified treatment of vector and tensor analysis shows the value of these 
subjects as tools in studying rigid-body dynamics, perfect fluids, differential 
geometry, finite deformation in elasticity, viscous fluids, compressible fluids, and 
turbulence. Excellent problems and exercises supplement the text. 


1955 . 341 pages . $7.00 


The Macmillan 


60 FIFTH AVENUE, NEW YORK 11, N.Y. 


Outstanding New Nc Graw-Hill Books 


NUMERICAL ANALYSIS Volume VI 


American Mathematical Society Proceedings of Symposia in Applied 
Mathematics 


Edited by Joun H. Curtiss. Ready for summer classes 


This new volume includes the papers presented at the symposium conference held under 
the above auspices at the Institute for Numerical Analysis at U.C.L.A., June 1953. The 
fundamental topics selected make this collection up to date and abreast with the increas- 
ing interest in and practical value of numerical analysis. The authors include Magnus R. 
Hestenes, Richard Bellman, and D. M. Lehmer. This volume will be of real value to 
applied mathematicians, computer-programming personnel, computer designers, and 
certain physicists and engineers. 


MODERN MATHEMATICS FOR THE ENGINEER 
By E. F. Becxenzacu, University of California, Los Angeles. Ready for fall classes 


The aim of this useful book is to generate in the minds of engineers and applied scientists 
engaged in research, design, and administration an awareness of the recent rapid ad- 
vancement in applied mathematical thought. This advancement resulted largely from 
the demands of modern engineering programming and design. It was made possible, in 
part, by recent advances in basic mathematics and statistics, and by the development of 
analog devices and digital computing machines of extremely high capacity and speed. 
Each chapter is written by an expert well known not only for his theoretical competence, 
but also for his applied experience. 


ADVANCED CALCULUS 


By R. Creicuton Buck, University of Wisconsin. /nternational Series in Pure and 

Applied Mathematics. Ready for Fall Classes 
An excellent graduate text for math majors, presupposing mastery of basic calculus and 
some differential equations. It gives a systematic and modern approach to the differ- 
ential and integral calculus of functions and transformations, and develops analytical 
techniques for attacking some of the typical problems which arise in applications of 
mathematics. Without retracing already familiar ground, it reviews elementary calculus 
with rigor. Includes 450 exercises. 


ENGINEERING ANALYSIS: A Survey of Numerical Procedures 


By StepHEN H. CranpALt, Massachusetts Institute of Technology. Ready for Fall 

Classes 
A second-year graduate text devoted to the advanced problems of engineering analysis. 
It considers briefly the formulation of mathematical models and the construction 
of computational programs for solving these with automatic or hand computing 
machines. The material is arranged according to a natural classification of the basic 
types of problems in engineering analysis, including equilibrium, eigenvalue, and propa- 
gation. It is the first book from an engineering viewpoint at this advanced level. 


McGraw-Hill Book Company, Inc. 


330 West 42nd Street New York 36, N.Y. 
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